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ABSTRACT 


Soil  erosion,  buffer  zone  deposition,  and  sediment  transport  processes  were 
studied  over  a  five-month  period  in  the  upper  reaches  of  Mill  Creek,  on  Fort  Jackson, 
South  Carolina.  The  objectives  of  the  study  were  to  approximate  the  erosional 
contribution  of  dirt  roads  to  sediment  detachment  and  transport,  to  assess  the  deposition 
of  eroded  soil  in  wetland  buffer  zones,  and  to  measure  suspended  sediment  leaving  the 
basin.  Repeat  cross-section  surveys  and  bulk  density  samples  of  selected  “rill  erosion” 
road  hillslopes  and  roadside  gullies  were  used  to  measure  soil  loss  during  the  period.  Soil 
loss  from  rainsplash  was  monitored  from  small  field  plots  (4.0  ft2).  Wetland  buffer  zone 
deposition  was  measured  by  a  series  of  plastic  mats  as  sediment  traps  across  a  buffer 
transect.  Particle-size  distributions  from  wetland  buffer  zone  deposits  were  measured  to 
evaluate  the  effectiveness  of  transport  and  storage  along  the  buffer  zone  transects. 
Suspended  sediment  concentrations  were  collected  from  two  tributaries  above  their 
confluence  and  from  the  main  channel  below  the  confluence. 

Combinations  of  instrumentation,  field  observations,  and  comparisons  with 
regional  values  reported  in  the  literature  and  from  model  output  were  used  to 
approximate  the  relative  importance  of  sediment  sources  and  sinks  in  the  upper  basin  of 
Mill  Creek.  Soil-loss  rates  and  sediment  deposition  results  were  used  to  evaluate  erosion, 
transport,  and  deposition  of  sediment  from  the  dirt  roads  and  buffer  zones. 

Unit  sediment  mass  from  dirt  road  small  box  plots  were  significantly  higher  than 
forested  small  box  plots  unit  sediment  masses.  Repeat  surveys  of  a  road  hillslope 
measured  a  continuous  degradation  of  the  slope  length  profile  and  a  net  soil  loss  of  15.68 
kg/m2.  Repeat  surveys  of  a  roadside  gully  measured  aggradation  and  degradation  of 


sediment  in  the  gully  and  a  net  soil  loss  of  1.42  kg/m2.  Sediment  deposition  was  found  to 
be  significantly  improved  with  buffer  zone  width.  During  moderate  storms,  deposition  of 
sediment  across  a  buffer  zone  transect  indicated  substantial  deposition  of  sediment  within 
the  first  ten  meters  of  the  buffer  zone  and  minimal  deposition  at  around  thirty  meters. 
Deposition  at  forty  meters  had  a  slight  increase  due  to  influence  from  overbank  flows 
from  streams. 

Suspended  sediment  concentrations,  stream  discharge,  and  sediment  flux  rates 
were  examined  (especially  sediment  -discharge  relationships)  to  identify  watershed 
responses.  Data  stratification  by  sub-basin,  by  rising  versus  falling  hydrograph  limbs,  by 
first  flush  phenomenon,  by  season,  and  by  land-use  changes  indicated  that  the 
relationship  between  suspended  sediment  and  discharge  at  the  main  channel  gage  site 
(Wl)  was  complex  due  to  varied  responses  from  the  two  tributaries.  A  first-flush  of 
sediment  was  dominant  in  the  smaller  northern  tributary  (WIN)  but  also  present  in  the 
larger  eastern  tributary  (W1E).  These  high  concentration  pulses  resulted  in  an  out-of¬ 
phase  relationship  between  peak  discharge  and  peak  suspended  sediment  concentrations 
causing  variance  in  rating  curves  and  complex  hysteresis.  Sediment  fluxes  from  five 
individual  storms  as  well  as  overall  unit  sediment  flux  and  unit  flow  rating  curves 
indicated  no  significant  difference  in  sediment  loadings  between  the  two  tributaries.  This 
is  attributed  to  the  importance  of  large  sediment  pulses  from  local  erosion-prone  sites  in 
both  tributaries.  In  short,  a  relatively  small  area  of  this  watershed  where  substantial 
sediment  contributions  are  generated  dominates  sediment  production  and  yield.  If  such 
sites  can  be  identified  and  controlled  NPS  pollution  could  be  greatly  reduced. 
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CHAPTER  I 


INTRODUCTION 

Concerns  with  non-point  source  (NPS)  pollution  have  been  fueled  by  recognition 
of  the  harmful  impacts  of  soil  loss,  channel  and  reservoir  in-filling,  degradation  of 
drinking  water  quality,  and  the  introduction  of  toxins  and  nutrients  to  water  sources. 
Realization  of  these  problems  has  made  the  identification  and  alleviation  of  such  sources 
imperative.  Section  319  of  the  Water  Quality  Act  of  1987  specifies  that  the  States  must 
implement  and  enforce  the  new  standards.  In  recent  years,  these  new  regulations  have 
increased  concern  for  and  research  of  NPS  pollution.  A  primary  type  of  NPS  pollution  is 
sediment.  Soil  erosion  not  only  removes  the  soil  surface  but  may  also  remove  much  of 
the  clay,  humus,  and  nutrients  from  the  soil  (James,  1998).  The  organic  layer  and  upper 
portion  of  the  mineral  horizons  of  soil  play  a  key  ecological  role.  These  are  the  layers 
from  which  vegetation  draws  much  nutrients  and  water.  High  rates  of  erosion  remove 
these  layers  faster  than  weathering  and  organic  decomposition  can  rebuild  them. 
Continued  removal  of  the  soil  A  horizon  reduces  the  infiltration  capacity  of  the  soil  and 
increases  the  amount  of  runoff,  which  further  increases  erosion  potential  (Dunne  and 
Leopold,  1978). 

A  study  area  was  selected  that  has  two  contrasting  tributaries,  one  tributary  that  is 
near  dirt  roads  and  firebreaks  and  a  second  tributary  is  located  in  the  center  of  a  wetland 
which  acts  as  a  buffer  zone  and  presumably  enhances  sediment  deposition  and  storage. 
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This  experimental  design  offers  an  opportunity  to  study  the  impact  of  roads  and  buffer 
zone  width  on  soil  erosion,  sediment  transport  and  deposition.  The  purpose  of  the  study 
was  to  assess  the  physical  properties  of  the  watershed  and  channel  system  affecting  the 
sediment  budget  of  the  Upper  Mill  Creek  watershed.  Limited  local-scale  erosion  data 
exist  for  southern  Sandhills  watersheds  especially  for  reforested  rural  watersheds  such  as 
upper  Mill  Creek  with  very  sandy  soils,  ample  relief,  and  abundant  dirt  roads.  This  study 
provides  quantitative  and  qualitative  insights  into  the  physical  processes  taking  place  and 
the  effectiveness  of  wetland  storage  and  dirt  road  erosion  in  the  Mill  Creek  watershed. 

Study  results  include  estimates  of  sediment  production,  deposition,  and  yield  from 
the  two  sub-basins.  Measurements  of  interrill  and  rill  erosion  rates  from  specific  dirt 
roads  and  forested  control  surfaces  are  used  to  determine  if  substantial  contributions  to 
the  overall  production  originate  from  the  dirt  roads.  Sediment  deposition  and  storage 
measurements  allow  the  identification  of  the  relative  importance  of  various  sediment 
storage  sites  and  overall  changes  in  storage  and  particle-size  distributions  in  this  type 
environment.  Sediment  yields  are  determined  at  the  base  of  the  two  tributaries  by 
sediment  sampling  and  stream  gage  recording.  While  an  ideal  study  would  include  a 
large  number  of  samples  of  each  of  these  erosion,  deposition,  and  yield  processes 
throughout  the  watershed,  limited  time  and  resources  available  for  this  study  severely 
constrained  the  number  of  samples  that  could  be  collected.  Yet,  it  was  thought  that  a 
proper  characterization  of  sediment  processes  in  the  basin  would  require  sampling  of  a 
diverse  array  of  sediment  source  and  storage  areas  as  well  as  transport  out  of  the 
watershed.  Thus,  this  study  was  designed  to  gather  a  limited  number  of  soil  and  sediment 
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samples  from  a  variety  of  environments,  to  characterize  the  watershed  processes  rather 
than  studying  a  single  factor  in  great  depth. 

SEDIMENT  PRODUCTION  AND  YIELD  PRINCIPLES 
Soil  Erosion 

Water  erosion  is  the  most  dominant  soil  erosion  process  in  the  southeastern 
United  States  and  includes  both  detachment  and  transport  of  sediment  (Colby,  1963). 

Lai,  1988,  defines  soil  erosion  as  “the  detachment  or  entrainment  of  soil  particles,  thus 
distinguishing  it  from  deposition  and  sedimentation  transport."  Sediment  transport  is  the 
actual  movement  of  sediment  for  example,  in  flowing  water  (Colby,  1963).  Sediment 
yield  is  the  removal  of  sediment  from  a  basin.  Water  moving  over  the  soil  is  a 
rudimentary  cause  of  soil  erosion.  To  generate  fluvial  erosion,  water  must  flow  over 
some  type  of  surface  that  allows  soil  particles  to  be  detached  and  begin  transport  down 
slope.  Important  factors  in  this  process  include  climate  (precipitation  and  vegetation), 
basin  size,  topography  (elevation  and  relief),  rock  type  (resistance  and  permeability  of  the 
surface),  and  human  activity  (land-use  and  management  pratices).  This  complex 
interaction  of  climate,  geology,  hydrology,  and  topography  effects  the  severity  of  erosion 
and  soil  loss  from  region  to  region  (USDA,  1983). 

Soil  erosion  is  a  natural  process,  but  if  the  magnitude  of  soil  erosion  is  accelerated 
then  harmful  effects  may  occur.  The  human  factor  in  soil  erosion  is  highly  variable 
depending  on  the  activity  and  the  environment  of  the  area.  Activities  such  as 
urbanization,  construction,  and  poor  agriculture  and  forestry  conservation  practices  may 
lead  to  severe  or  accelerated  erosion  (Ritter,  Kochel,  and  Miller,  1995). 
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Soil  erodibility  is  a  measure  of  the  vulnerability  of  soil  to  detachment  and 
movement  (Lai,  1988).  Rock  type  characteristics  effect  the  erodibility  and  infiltration 
capacities  that  are  responsible  for  sediment  production  and  runoff  amounts  of  a  basin. 

Soil  characteristics  such  as  particle-size  distribution,  percentage  organic  content,  soil 
structure  and  soil  permeability  are  important  factors  in  soil  erodibility  (Weaver  and 
Lineback,  1981). 

Hillside  Erosion  Processes 
Rainsplash  Erosion 

As  rainfall  strikes  the  surface,  it  contains  a  relatively  large  amount  of  kinetic 
energy.  The  applied  force  by  raindrops  acts  upon  the  soil  surface  to  detach  sediment 
particles  and  begin  the  transport  downslope.  This  initial  detachment  is  “ rainsplash 
erosion ”  (Troeh  et  al.,  1991)  which  depends  on  numerous  soil  and  basin  characteristics. 
The  most  notable  characteristic  is  the  soil  erodibility,  which  results  from  many  inherent 
properties  and  changes  as  the  soil,  reacts  to  climate,  biota,  and  land-use  changes.  Soil 
texture  influences  the  detachment  and  transport  of  sediment  particles.  Course  sand 
particles  resist  transport  while  silty  and  clayey  soils  have  strong  cohesive  properties  that 
help  resist  detachment.  As  the  finer  grain  soils  detach  and  deposit,  they  may  form  a  crust 
or  impermeable  layer  which  may  increase  surface  runoff  (Lai,  1988) 

Sheet  Wash  Erosion 

As  rainfall  intensity  exceeds  infiltration  capacity  or  soil  becomes  saturated,  and 
water  begins  to  run  off,  a  different  type  of  erosion  takes  place.  Overland  flow  erodes  the 
soil  surface  over  time  by  “sheet  wash ”  or  “ interrill  erosion ”  (Dunne  and  Leopold,  1978). 
“ Rainsplash  erosion ”  influences  the  amount  of  sheet  erosion  by  entraining  particles.  The 
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physical  properties  of  particle  sizes,  cohesiveness,  porosity,  antecedent  moisture, 
vegetation,  slope  gradient,  and  slope  length  effect  sediment  production  and  the  ability  of 
runoff  to  detach  and  transport  sediment  (Colby,  1963).  If  uniform  removal  of  the  soil 
surface  by  runoff  does  not  occur,  small  channels  driven  by  shear  stresses  of  channel  flows 
may  form  in  areas  of  least  resistance.  This  process  of  “ rill  erosion ”,  causes  an  increase  in 
the  efficiency  and  intensity  of  erosion  and  transport  (Dunne  and  Leopold,  1978).  These 
two  different  types  of  erosion  play  an  important  role  in  modeling  soil  erosion.  Each 
process  is  dominated  by  different  forces  and  requires  different  calculation  methods.  Each 
process  also  leaves  different  evidence  in  the  field. 

The  majority  of  the  sediment  generated  on  hillslopes  stays  within  the  basin  due  to 
deposition  and  storage  (ASCE,  1975;  Lai,  1994;  Meade,  1982).  The  topography  and 
vegetation  cover  of  the  area  greatly  influences  this  process.  Concave  slopes,  crop  field 
boundaries,  floodplains,  and  reservoirs  act  as  collection  points  for  sediments  (Lai,  1988). 
Sediment  yields  from  a  watershed  often  account  for  only  a  small  portion  of  sediment  loss 
by  erosion.  These  deposition  areas  reduce  the  energy  of  the  channeled  flow  that  is 
transporting  the  sediment  particles.  Decreased  gradient,  sloe-length,  changes  in  channel 
geometry,  and  increased  friction  due  to  rougher  terrain  and  vegetation  all  reduce  the  flow 
velocity.  This  reduction  in  velocity  allows  larger  particles  to  sink  due  to  gravity 
(Bagnold,  1973). 

Evidence  of  the  role  of  basin  size  in  soil  loss  and  sediment  transport  indicates  that 
the  sediment  yield  from  a  watershed  is  greatest  in  smaller  basins.  These  generalizations 
suggest  that  most  sediment  is  produced  in  the  upper  reaches  and  a  major  portion  of  it  is 
stored  in  the  floodplain  to  be  removed  through  geomorphologic  and  hydrologic  processes 
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(Roehl,  1962).  The  “ sediment  delivery  ratio ”  of  a  watershed;  that  is,  the  ratio  of  sediment 
yield  to  sediment  production,  is  an  assessment  tool  of  land-use  practices  and  often  an 
objective  of  erosion  models  (Dunne  and  Leopold,  1978;  Lai,  1994)  (Figure  1-1). 
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Figure  1-1.  Example,  SDR  verses  Drainage  Area 
Source:  USDA,  SCS  National  Engineering  Handbook,  1975 
Deposition  in  Buffer  Zones 

Buffer  zones  are  natural  filters  and  traps  for  both  nutrients  and  sediment  that 
impose  important  controls  for  improving  water  quality,  and  have  been  shown  to  be 
valuable  in  reducing  the  amount  of  (NPS)  pollution  to  receiving  streams.  A  United  States 
Ddepartment  of  Agriculture  (USDA)  recommended  average  buffer  zone  width  is  50  m 
(Lowrance  et  al.,  1997).  One  type  of  buffer  zone  is  a  riparian  buffer.  The  USDA  defines 
this  as  a  belt  of  native  trees  and  shrubs  located  adjacent  to  and  upslope  from  bodies  of 
water.  Riparian  zones  consist  of  three  zones.  Zone  one  is  typically  5  to  7  meters  (15-20 
feet)  and  consists  of  streamside  vegetation.  Zone  two  begins  at  the  end  of  zone  one  and 
extends  another  7  to  20  meters  (20-60  feet).  It  is  a  transitional  forest  that  includes 
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conifers,  shrubs,  and  some  hardwoods.  Timbering  may  be  done  in  zone  two  if  soil 
stability  is  not  adversely  effected  and  zone  one  is  not  compromised.  Ruffin  (1998) 
recommends  that  only  50%  of  the  timber  should  be  removed.  Zone  three  is  used  if  the 
site  is  next  to  tilled  or  grazed  land  and  it  should  extend  at  least  7  meters  (20  feet)  past 
zone  two.  A  general  rule  of  thumb  is  that  the  width  of  zones  one  and  two  combined 
should  be  one-third  the  distance  to  the  farthest  sediment  source  area  (Weik,  1999). 
WEPP  Hillslope  Erosion  Model 

Through  past  examination  of  standard  plots  and  applying  the  related  plot  data  to 
small  watersheds  and  field  sites,  erosion  prediction  models  such  as  the  Water  Erosion 
Prediction  Project  (WEPP)  is  a  tool  for  researching  the  physical  processes  involved  in 
sediment  erosion,  storage,  and  transport.  By  utilizing  climate,  topography,  soil 
characteristics,  vegetation,  and  human  activities,  careful  application  of  the  WEPP  model 
may  assist  in  land-use  management. 

The  WEPP  is  a  physically-based  model  that  predicts  erosion,  downslope 
deposition,  and  sediment  yield  on  a  daily  time  step.  It  can  estimate  infiltration,  interrill 
erosion,  and  rill  erosion  to  predict  runoff  and  sediment  yield  from  an  elevated  area 
downslope  at  a  local  scale  (Morfin  et  al.,  1996).  The  dominant  forces  are  interrill 
(raindrop  splash  and  sheet  flow)  erosion  and  rill  erosion  (channelized  flow).  Effective 
hydraulic  conductivity  (Ke)  is  a  major  factor  in  WEPP  for  calcualting  runoff  (Elliot, 
Foltz,  andRembolt,  1994). 

WEPP  has  several  conceptual  parameters  that  estimate  soil  detachment  and 
deposition.  The  four  main  input  files  are  management,  soil,  slope,  and  climate  (Morfin, 
et  ah,  1996).  The  management  input  file  allows  for  descriptions  of  vegetation,  human 
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management  practices,  and  initial  conditions  prior  to  an  event  or  time  simulation.  The 
soil-input  file  describes  soil  texture,  albedo,  percentage  saturation,  interill  erodibility  (Ki), 
rill  erodibility  (Kr),  and  critical  shear  of  the  soil  (tc).  It  can  have  up  to  ten  soil  map  layers 
(along  the  hillslope  profile)  and  can  extend  up  to  2  meters  in  depth.  The  slope-input  file 
specifies  the  length  and  width  of  hillslopes.  A  slope  length  can  be  modeled  in  ten 
different  gradients  enabling  a  detailed  description  of  complex  topography.  The  climate- 
input  file  describes  the  daily  maximum,  minimum,  and  dew  point  temperatures,  rainfall 
intensity  and  duration,  and  discharge  time  to  peak,  solar  radiation,  and  wind  speed  and 
direction  (Elliot,  Foltz,  and  Rembolt,  1994). 

STUDY  AREA  DESCRIPTION 

Upper  Mill  Creek  is  located  east  of  Columbia,  South  Carolina  on  the  Fort  Jackson 
Military  Reservation  (Figure  1-2).  Columbia  is  located  in  the  Sandhills  between  the 
Piedmont  and  the  Upper  Coastal  Plain.  The  climate  in  the  Midlands  is  relatively  mild 
with  average  annual  temperatures  of  approximately  14  °C  (low  60s  °F)  and  an  average 
annual  precipitation  ranging  from  106  to  120  cm/yr  (42-47  inches/yr).  From  February 
1999  to  January  2000  the  study  area  suffered  from  an  extreme  drought  resulting  in  record 
low  stream  flows  (SC  DNR  unpublished  data,  2000).  For  this  study,  storm  events 
occurring  between  April  through  September  will  be  categorized  as  summer  storms  and 
the  storm  events  occurring  between  October  through  March  will  be  categorized  as  winter 
storms. 
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Figure  1-2.  Upper  Mill  Creek  Study  Area 


The  total  area  of  the  watershed  is  6.6  km2.  The  study  watershed  is  in  the  upper 
reaches  of  the  west  tributary  of  Mill  Creek  (henceforth  referred  to  as  upper  Mill  Creek  or 
Wl).  W1  is  mainly  forested,  with  82.7%  woodland  area,  1.8%  grassland  area,  11.7% 
wetlands,  and  3.8%  dirt  roads  (Dean  et  al.,  1998).  An  extensive  system  of  dirt  roads  and 
firebreaks  borders  and  crosses  the  entire  drainage  area  (Figure  1-3).  Based  on  a 
preliminary  examination  of  historical  aerial  photographs,  these  fire  breaks,  dirt  roads,  and 
range  complexes  were  first  introduced  between  1955  and  1959.  This  dense  network  of 
dirt  roads  and  firebreaks  appears  to  accelerate  soil  erosion  and  act  a  source  of  sediment. 
The  channel  network  largely  reflects  roadside  ditches  mapped  in  the  field  in  1998  (James, 
personal  communication)  and  updated  by  this  study.  (Figures  1-2  and  1-3). 


9 


Figure  1-3.  Study  Area  Dirt  Road  and  Firebreak  Network 
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Figure  1-4.  WIN  (left)  and  W1E  (right)  Confluence  at  W1  Gage  Site 

(View  Upstream) 

The  main  channel  bifurcates  a  few  meters  above  the  W 1  outlet  into  two  main 
tributaries  (WIN  and  W1E)  (Figure  1-4).  Both  channels  are  ditched  and  straight  in  their 
lower  reaches  (near  the  W1  gage).  The  north  tributary  (WIN)  is  paralleled  by  a  dirt  road 
and  has  direct  road  runoff  entering  WIN.  WIN  has  a  substantially  narrower  buffer  zone 
along  the  channel  compared  toWlE.  The  W IE  channel  ends  about  100  m  upstream  of 
the  outlet  within  a  broad  wetland  buffer  zone  approximately  three  times  as  wide  as  WIN. 
W1E  has  a  large  tributary  that  drains  from  an  active  rifle  range  that  supplies  a  large 
volume  of  runoff  and  sediment.  This  tributary  at  one  location  flows  directly  on  a  dirt 
road  for  100  m  just  prior  to  joining  the  W IE  channel  (Figure  1-5). 
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Figure  1-5.  W1E  Channel  Diversion  Flowing  on  a  Dirt  Road 

(View  Upstream) 

WIN  comprises  about  20%  of  the  total  drainage  area  of  the  W1  watershed  (Figure 
1-6).  It  has  an  earthen  dry  dam  about  1000  m  upstream  from  the  W1  outlet.  The  earthen 
dam  has  no  spillway  but  a  24  in.  culvert  allows  sediment  to  be  discharged  through  the 
reservoir.  The  WIN  ditched  channel  bottom  between  the  dirt  roads  and  W1  oulet  has 
shown  no  sign  of  aggradation  or  degradation  during  the  study  period.  WIN  crosses  a  dirt 
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road  and  flows  through  a  gently  sloping  transition  into  the  wetland.  From  approximately 
100  m  into  the  wetland  WIN  is  ditched  to  its  confluence  with  W1E  at  the  Wloulet. 
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Figure  1-6.  WIN  and  W1E  Drainage  Basins  within  the  Upper-most  West 

Tributary  of  Mill  Creek 

W1E  comprises  about  80%  of  the  total  drainage  area  of  the  W1  watershed  (Figure 
1-6).  W IE  appears  to  be  stable  near  the  confluence,  although  some  upstream  reaches  in 
the  wetland  have  aggraded.  The  head  of  the  W1E  channel  does  not  appear  to  be  eroding 
but  unlike  WIN,  comes  to  an  abrupt  end  about  30  m  below  a  dirt  road  named  Old 
Hartsville  Guard  Road. 
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Extensive  timber  harvesting  occurred  in  early  December  and  other  assoicated 
changes  in  the  W1  watershed  included  road  grading,  increased  vehicular  activity,  the 
thinning  of  pine  trees,  and  the  addition  of  timber  debris  unto  roads  and  in  the  WIN  and 
W1E  channels. 

Previous  Studies 

Studies  conducted  in  1997  and  1998  (Dean  et  al.,  1998;  2000)  provided  a  base 
line  of  stream  flow  and  suspended  sediment  discharge  data  for  the  main  W1  channel 
below  the  confluence  of  W1E  and  WIN.  These  studies  indicated  that  for  high-intensity 
storms  with  high  antecedent  soil  moisture,  suspended  sediment  concentrations  increased 
quickly  in  response  to  storm  events.  During  an  intense  convective  thunderstorm  on  29 
July  1997,  suspended  sediment  concentrations  peaked  early  in  the  event  displaying  the 
"first  flush  phenomenon"  (from  the  WIN  channel).  A  second  sediment  concentration 
peak  was  observed  approximately  2.25  hours  later  and  corresponded  to  the  time  to  peak 
concentration  of  the  W1E  channel.  A  less  intense  frontal  rainfall  event  on  26  Oct  1997 
produced  different  results.  This  storm  did  not  produce  the  "first  flush  phenomenon"  seen 
in  the  high  intensity  storm  but  had  a  slower  time  of  peak  concentration  with  lower 
sediment  fluxes.  There  was  also  a  second  peak  sediment  concentration  observed  for  this 
storm  (Dean  et  al.,  1998;  Atkins,  personal  communication,  1999).  That  study  suggested 
that  storm  type,  antecedent  conditions,  and  the  physical  nature  of  the  basin  effects  the 
WIN  and  W1E  suspended  sediment  concentrations.  The  fast  response  with  high 
concentrations  from  W  IN  suggests  inadequate  buffer  zones  and  quick  transport  of 
sediment  particles  along  the  dirt  roads.  The  slower  response  from  the  east  channel 
indicates  better  buffering,  longer  travel  distances,  and  slower  time  to  conentrations. 
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OBJECTIVES 


The  W 1  watershed  provides  a  good  study  area  to  observe  the  physical  processes 
involved  in  hillslope  soil  erosion  and  sediment  transportation  in  a  Sandhills  environment. 
The  bifurcation  of  W 1  into  two  distinct  sub-basins  enables  a  study  design  to  compare  the 
effects  of  soils,  vegetation,  buffer  zone  characteristics  and  land  disturbing  activities 
(roads,  firebreaks,  and  timber  harvesting)  on  the  erosion  and  delivery  of  sediment. 

Several  hypotheses  can  be  generated  and  tested  concerning  these  different 
responses  or  effects.  Hypothesis  1  is  that  most  entrainment  and  transport  of  sediment 
occurs  along  the  dirt  roads  as  compared  to  forested  areas.  Hypothesis  2  is  that  the 
wetland  buffer  stores  a  substantial  portion  of  the  total  sediment  generated  resulting  in  a 
reduction  of  sediment  transport  and  delivery  to  the  stream.  Hypothesis  3  is  that  WIN  will 
have  a  higher  sediment  yield  per  unit  area  than  W1E,  because  (A)  it  has  a  much  smaller 
buffer  zone  along  its  channel  compared  to  W1E  and  because  (B)  the  proximity  of  the  road 
to  the  channel  enhances  sediment  delivery  (Phillips,  1989). 
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CHAPTER  II 


FIELD,  LABORATORY,  AND  STATISTICAL  METHODS 

A  wide  variety  of  methods  were  used  in  this  thesis  in  order  to  evaluate  a  range  of 
processes  and  environments.  This  chapter  which  is  divided  into  three  sections:  field, 
laboratory,  and  statistical  methods.  Field  methods  include  stream  discharge 
measurements,  suspended  sediment  sampling,  small  box  plot  interrill  erosion 
measurements,  repeat  road  and  gully  cross-sectional  surveys  for  rill  erosion,  bulk  density 
samples,  and  sediment  deposition  measurements  using  artificial  grass-mats.  Laboratory 
methods  utilized  were  suspended  sediment  concentrations  through  filtration  bulk  density, 
removal  of  organic  matter  by  loss-on-ignition  (LOI),  sediment  textural  analysis  through 
wet  and  sonic  sieving.  Statistical  methods  include  regression,  t-test,  and  ANCOYA. 
FIELD  METHODS 

Several  erosion  processes  were  monitored  through  field  observation  and 
measurements.  Discussion  of  field  methods  will  begin  with  rainfall,  then  progress  to 
runoff,  erosion,  sediment  deposition,  and  sediment  yield. 

Rainfall 

Rainfall  data  were  recorded  at  the  Bravo  9  tipping  bucket  gauge  site, 
approximately  1.5  kilometers  southeast  of  the  study  watershed’s  outlet  (G.  Carbone, 
unpublished  data).  Total  rainfall  for  the  study  was  61.5  cm  with  an  annual  total  of  74  cm, 
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most  of  which  was  observed  during  the  period  from  late  June  to  September  1999. 


9B  Precip  1999 


Figure  2-1.  (1999  Daily  Precipitation  at  Bravo  9  Raingage  (G.  Carbone, 

Unpublished  data)). 

Runoff 

Previous  studies  at  W1  established  discharge  functions  using  data  from 
continuous  stage  recorders  and  streamflow  discharge  measurements  of  discharge  to 
establish  stage-discharge  relationships  (Dean  et  al  1998,  A.  James  personal 
communication,  1998).  Sensitivity  of  the  stage  recorders  was  within  +  1  cm  following 
calibration.  The  existing  W1  stage  sensor  and  data  logger  has  had  a  continuous  record 
since  April  1997.  To  measure  stages  at  the  W1E  channel  outlet  a  stage  sensor  and  data 
loggers  were  installed  on  May  19, 1999  to  record  at  fifteen-minute  intervals.  The  initial 
stage  sensor  failed  and  was  reestablished  June  1999.  This  instrumentation  facilitated  the 
calculation  of  runoff  from  stage  discharge  relationships  developed  by  discharge 
measurements.  Staff  gages  were  attached  by  rigid  wire  to  metal  stakes  secured  in  the 
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channel  bottom.  Stage  recorders  were  anchored  with  rigid  wire  to  plastic  pipe  stilling 
wells  in  the  W1E  tributary  and  the  W1  main  channel  and  anchored  to  metal  stakes. 


Table  2-2.  Sketch  of  Gage  Site  at  W1  Outlet  and  WIN  and  W1E  Confluence. 


Discharge  measurements  were  taken  from  a  cross-section  at  the  W1E  sampling 
site  and  related  to  stages  at  staff  gage  7  (G7)  at  this  site.  The  cross-section  was  divided 
into  equal  5  cm  horizontal  segments  using  a  board  that  was  staked  across  the  channel. 
There  were  at  least  20  sections  in  all  discharge  measurements.  Given  the  symmetrical 
shape  and  uniform  depth  of  the  channel  discharge  is  fairly  uniformly  distributed  across 
the  section.  Thus,  no  single  measurement  composed  a  large  percentage  of  the  total.  A 
depth  reading  was  taken  with  a  standard  wading  rod  utilizing  the  .6  method  (Dunne  and 
Leopold,  1978).  Velocities  were  measured  with  a  Marsh-McBimey  digital  flow  meter. 
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The  segment  areas  were  calculated  with  the  velocity  measurement  in  the  segment  center. 
The  segment  discharges  were  calculated  and  summed  to  get  the  channel  discharge  (e.g. 
Table  2-1).  The  stage  at  G7  was  recorded  during  discharge  measurement  to  establish  a 
stage-discharge  rating  curve  for  the  W1E  tributary  (Table  2-2). 

In  order  to  convert  large  amounts  of  stage  data  from  the  data  loggers  to 
corresponding  discharge  the  data  in  Table  2-2  were  functionalized  by  regressing  the  G7 
gage  height  against  discharge  using  a  second-order  polynomial  (Figure  2-3).  This 
relationship  is  strong  (R2  =  0.98)  within  the  limited  range  of  flows  measured.  Higher 
discharges  need  to  be  measured  to  validate  the  relationship  in  Figure  2-3.  Due  to 
conditions  during  the  study  period  there  were  limited  opportunities  to  measure  high  flows 
and  other  field  activities  often  competed  for  attention  during  those  events.  Thus,  the 
discharge  for  W1  was  calculated  using  the  stage-discharge  regression  from  a  previous 
study  (Figure  2-4)  (Dean  et  al.,  2000). 

Given  that  the  W1  stream  gage  and  sediment  sampling  site  is  about  three  meters 
below  the  confluence  of  the  two  tributaries,  it  is  assumed  that  changes  in  storage  between 
the  two  gages  are  negligible.  Therefore  discharge  rates  of  water  and  sediment  passing 
through  the  two  tributaries  should  sum  to  the  totals  observed  as  below  the  confluence. 
Thus  a  simple  relationship  follows  in  relation  to  discharge  at  Wl: 

Qwi  =  Qwin  +  Qwie  (Equation  1) 

Where,  Q  is  discharge  (1  /  s)  and  subscripts  Wl,  WIN,  and  W1E  refer  to  gage  sites  at  the 
Wl  outlet,  north  tributary,  and  east  tributary,  respectively.  These  assumptions  allowed 
calculation  of  Qwin  by  subtracting  Qwie  from  Qwi  instantaneous  readings  taken  within  a 
three-minute  sampling  period.  Changes  in  flow  are  typically  minimal  at  this  time  scale  so 


19 


Qwin  estimates  should  be  accurate  within  the  precision  limits  of  the  individual  discharge 


measurements. 


W1EQ  ~  | 

26-Jun-99 

interval  (Measured  Right  to  Left-Looking  upstream)  1 

Interval  (m) 

Depth  (m) 

Sect  Veiocity  (m/s) 

Sect  Width  (m) 

Sect  Q  (m3/s) 

0.25 

0.00 

0.00 

0.30 

0.14 

-0.01 

0.025 

0.0035 

-0.00002 

0.35 

0.14 

-0.01 

0.05 

0.0070 

-0.00006 

0.40 

0.20 

0.00 

0.05 

0.0100 

-0.00003 

0.45 

0.22 

0.01 

0.05 

0.0110 

0.00007 

0.50 

0.26 

0.01 

0.05 

0.0130 

0.00012 

0.55 

0.27 

0.01 

0.05 

0.0135 

0.00012 

0.60 

0.02 

0.05 

0.0145 

0.00027 

0.65 

0.29 

0.02 

0.05 

0.0145 

0.00035 

0.70 

0.29 

0.05 

0.05 

0.0145 

0.00071 

0.75 

0.32 

0.06 

0.05 

0.0160 

0.00102 

0.80 

0.33 

0.07 

0.05 

0.0165 

0.00116 

0.85 

0.33 

0.06 

0.05 

0.0165 

0.00106 

0.90 

6.34 

0.07 

0.05 

0.0170 

0.00119 

0.95 

0.34 

0.07 

0.05 

6.0170 

0.00114 

1.00 

0.33 

0.06 

0.05 

0.0165 

0.00101 

1.05 

0.32 

0.04 

0.05 

0.0160 

0.00068 

1.10 

0.30 

0.03 

0.05 

0.0150 

0.00050 

1.15 

0.29 

0.02 

0.05 

0.0145 

0.00035 

1.20 

0.27 

0.02 

0.05 

0.0135 

0.00033 

"T25 

0.24 

0.00 

0.025 

0.0060 

0.00002 

1.30 

0.24 

0.00 

0.00998 

Stage  Gage  7  (mm): 

310 

Total  Q  (L/s)= 

r“  9-98 

Table  2-1.  Example  W1E  Discharge  Calculation 


Date 

G7(MM) 

Q(l/s) 

15-Jun 

270 

5.83 

16-Jun 

354 

21.94 

26-Jun 

310 

9.98 

9-Sep 

260 

3.63 

10-Sep 

256 

3.49 

15-Sep 

284 

8.82 

15-Sep 

310 

11.32 

15-Sep 

305 

8.45 

28-Sep 

375 

25.45 

28-Sep 

355 

21.05 

Table  2-2.  Stage  Discharge  readings  at  W1E 
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Q  (I/s) 


Figure  2-4.  W1  Stage-Discharge  Rating  Curve  (Dean  et  al.,  1998) 
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Sediment  Production 


Interrill  box  locations  were  constructed  at  four  sites:  Rl,  near  the  outlet  on  the 
WIN  tributary,  R2  upstream  of  the  W IE  tributary  in  proximity  to  the  clay  pit,  R3 
upstream  on  the  WIN  tributary  near  the  check-dam,  and  R4  upstream  of  the  W1E 
tributary  near  the  W1E  large  tributary  (Figure  2-5).  Merrill  sediment  was  measured 
using  a  small-field  box  plots  described  in  Lai  (1994).  Each  Merrill  field  plot  consisted  of 
a  4.0  ft2  wooden  box  which  was  placed  on  sites  representative  of  the  dirt  roads  or  in  the 
forest  as  controls.  The  box  has  a  spout  on  its  lower  end  that  emptied  runoff  and  sediment 
into  a  500  ml  collector  (Figure  2-6). 


Figure  2-5.  Merrill  and  Rill  Erosion  Sites 
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Figure  2-6.  Small  Box  Plot  (Merrill  Box) 

Rill  erosion  sites  were  at  the  R3  and  R2  sites,  the  hillslope  profile  was  at  R3,  and 
the  gully  cross-section  was  at  R2  (Figure  2-5).  Erosion  and  deposition  along  dirt  road 
rills  was  measured  with  cross-section  profiles  across  rills  at  monumented  locations  (e.g. 
Sirvent  et  al.,  1995).  On  the  dirt  road  hillslope  (R3  site)  (Figure  2-7)  four  cross-sections 
were  measured  (0  m,  15  m,  45  m,  and  75  m)  along  with  a  slope  profile  down  the  middle 
of  the  road  at  0  m,  15  m,  35  m,  60  m,  75  m,  and  98  m.  The  cross-sections  were  surveyed 


with  a  manual  transect  using,  100  m  surveying  tapes,  and  a  7  m  survey  pole  to  take 
readings  at  about  5-cm  intervals.  The  R2  site  was  a  roadside  gully  approximately  1  m 
wide  with  a  slope  of  2  percent  (Figure  2-8).  Two  cross-sections  approximately  2.3  m 
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across  and  2.6  m  apart  were  measured  at  1-cm  intervals  by  placing  a  level  line  across  the 
tops  of  monumented  survey  stakes  and  measuring  the  cross-section  depth  with  a  metric 
ruler.  Cross-section  areas  were  calculated  using  the  same  procedure  as  in  the  discharge 
procedure.  Erosion  volumes  were  calculated  as  the  product  of  the  difference  in  areas  and 
the  hillslope  distance  between  cross-sections. 

Bulk  density  samples  were  collected  using  a  100-cm3  cylinder.  Bulk  density  was 
used  to  convert  sediment  volumes  to  mass  to  assess  the  amount  of  sediment  being 
entrained  and  deposited.  Samples  were  taken  along  the  soil  surface  in  areas  with  no 
pebbles  or  macro  organics  at  approximately  5  cm  depths.  Two  samples  were  taken  in  the 
side  spoil  generated  from  recent  grading  and  three  samples  were  taken  from  within  rills 
representative  of  surface  along  the  hillslope. 


Figure  2-8.  R2  Gully  Survey  Location.  Dirt  Road  north  of  Old  Hartsville 
Guard  Road  that  empties  into  W1E  (December  1999) 

Deposition 

Wetland  buffer-zone  deposition  was  measured  utilizing  plastic  mats  (Figure  2-9) 
systematically  placed  in  representative  areas  of  likely  deposition  below  roads  at  wetland 
boundaries.  Depositional  mats  were  placed  in  the  same  general  location  as  the  four 
erosion  sites.  Ml  was  along  WIN  above  the  ditched  channel  at  the  edge  of  the  wetland, 
M2  was  above  W1E  below  the  sand  pit,  M3  was  on  WIN  below  the  earthen  check  dam, 
and  M4  was  above  the  W1E  ditched  channel  (Figure  2-10).  Mats  were  placed  at  intervals 


between  10  m  -  15  m  inside  the  buffer  zone  starting  from  the  road  and  toward  the 
channel  and  secured  with  six-inch  turf  spikes.  Depositional  mat  transects  topography  had 
only  slight  relief  with  slopes  less  than  two  percent.  Vegetation  varied  little  between  sites 
and  consisted  of  mainly  pines,  hardwoods,  evergreen  shrubs,  and  grasses.  Each  mat  had 
an  area  of  0.297  m2  and  was  a  rigid  plastic  turf  that  had  a  roughness  value  higher  than  the 
surrounding  surfaces  due  to  the  vertical  rigidity  and  density  of  the  1-  cm  artificial  grass 
blades.  Additionally,  the  secured  mat  was  not  washed  aside  like  the  pinestraw,  leaves, 
and  grasses  that  dominated  the  litter  layer.  Thus  sediment  deposits  may  have  been  greater 
than  on  adjoining  areas.  Nevertheless,  spatial  patterns  of  deposits  should  be 
representative  of  natural  surfaces.  Field  observations  after  storm  events  found  little 
visual  difference  between  texture  of  sediment  deposits  around  the  mats  trapped  by 
underlying  grasses,  roots,  and  sticks  and  sediment  deposited  on  the  mat. 


Figure  2-10.  Depositional  Mat  Photo.  Ruler  is  one  Foot  in  Length 
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Deposition  Mat  Study  Sites 


Figure  2-10.  Depositional  Mat  Locations 


27 


Sediment  yields 


Manual  depth  integrated  (grab)  methods  (Hadley  and  Walling,  1984)  were  used  to 
collect  500  ml  suspended  sediment  samples,  at  WIN,  W1E,  and  Wl.  In  addition,  an 
automated  ISCO  sampler  (ISCO,  1996)  was  used  to  sampled  at  Wl  with  the  intake 
positioned  in  the  center  of  the  channel  for  the  intermittent  periods  through  August  1999. 
These  techniques  allow  for  representative  samples  of  suspended  sediment  particles 
(Horowitz,  1991).  Manual  sampling  followed  guidelines  in  the  USDA,  SCS,  National 
Engineering  Handbook,  Chapter  4, 1975. 

Given  the  short  distances  between  gages  (Figure  1-4)  and  the  small  amount  of  bed 
material  in  the  channels,  changes  in  sediment  storage  between  the  Wl,  W1E,  and  WIN 
sediment  sampling  sites  were  assumed  negligible.  Therefore,  sediment  yields  at  the 
mouth  of  the  two  tributaries  were  assumed  to  sum  to  the  total  yield  below  the  confluence: 

Qswi  =  Qswin  +  Qswie  (Equation  2) 

Where  Qs  is  sediment  discharge  (g  /  s).  This  is  similar  to  the  assumption  made  for  runoff 
(Equation  1)  except  that  sediment  tends  to  travel  in  pulses  so  the  three-minute  time 
difference  between  the  samples  could  be  more  critical.  To  test  validity  of  the  relationship 
assumed  by  equation  2  sediment  samples  were  collected  at  all  three  gages  and  discharges 
calculated  at  all  three  sites. 

LABORATORY  METHODS 
Bulk  Density 

Bulk  density  is  the  ratio  of  dried  mass  to  volume  (g/cm  ).  Bulk  density  was 
determined  for  erosion  cross-section  samples  by  drying  the  100  cc  samples  overnight  at 
103  °C,  cooling  in  a  dessicator,  and  weighing  (Blake  and  Hartge,  1986).  Collected 
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samples  had  little  to  no  macro  organics,  which  were,  therefore,  included  in  the  overall 
weight. 

Deposition  Mat  Sediment 

Mats  were  carefully  lifted  in  the  field  and  transported  to  the  lab  in  plastic  bags. 
Larger  particles  were  removed  to  a  500-ml  beaker  and  bags  and  mats  were  thoroughly 
rinsed  into  a  5-gallon  bucket.  The  sample  was  allowed  to  settle  for  72  hours,  the  excess 
water  was  drained,  the  sediment  moved  to  a  labeled  beaker,  dried  at  103  °C  for  24  hours, 
and  weighed  using  an  analytical  balance.  A  unit  sediment  mass  was  obtained  by  dividing 
the  total  dried  weight  by  the  mat  area.  The  sediment  was  then  stored  in  a  dessicator  prior 
to  particle  size  analysis. 

Organic  Matter 

To  determine  the  mass  of  organic  matter  in  the  mat-deposition  samples,  the  loss- 
on-ignition  (LOI)  method  was  used  as  described  in  the  Standard  Methods  for  the 
Examination  of  Water  and  Wastewater  (Eaton  et  al.,  1995).  Samples  were  placed  in  a 
muffle  furnace  and  heated  at  550  °C  for  forty  minutes.  The  sample  was  removed,  allowed 
to  cool  in  a  dessicator,  and  weighed.  The  percentage  organic  matter  lost  was  calculated 
as: 

%OM  =  Wh-Wi  *  100  (Equation  3) 

Wd 
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where  %  OM  is  percent  organic  matter,  Wd  is  mass  of  dried  sample,  and  Wi  is  mass  of  the 
ignited  sample. 

Particle-Size  Distribution  on  Deposition  Mats 

Wet-sieving 

Once  the  sediments  from  the  deposition  mats  had  been  dried  and  the  organic 

matter  removed  by  LOI,  the  samples  were  wet-sieved  and  sonic  sifted  to  determine  the 

fine  and  sand  fractions.  The  sample  was  carefully  ground  in  a  mortar  to  break  up 

aggregates  using  a  rubber-tipped  pestle.  From  the  sample,  30  grams  (or  the  entire  sample 

if  less  than  30  grams)  was  remove  with  a  plastic  spoon  and  placed  in  a  63um  sieve,  and 

fines  were  washed  out  by  a  steady  stream  of  water.  The  remaining  material  was 

extracted,  placed  in  a  beaker,  oven  dried  at  103°C  for  24  hours,  cooled,  and  weighed  to 

get  the  total  mass  of  sand.  The  difference  from  the  original  mass  was  the  percent  fines: 

%F  =  Wh-Ws  *  100  (Equation  4) 

Wd 

where  %  F  is  percent  fines  and  Ws  is  weight  of  sieved  sediment.  The  remaining  material 
was  then  sonic  sifted  to  get  the  sand  distribution. 

Sonic  Sifting 

An  ATM  sonic-sifter  was  used  to  analysis  sand  particle-size  distributions  (ATM, 
no  date).  The  sieves  were  cleaned  by  gently  tapping  each  individual  sieve.  Six  plastic 
sieves  with  a  wire  mess  of  decreasing  geometric  size  from  2000  um  to  63  um  were  used. 
The  sieves  were  stacked  in  the  following  order:  2000, 1000,  500,  250, 125,  and  63  um. 
Approximately  15  -20  grams  of  dry  wet  sieved  material  was  placed  in  the  top  sieve  and 
inserted  into  the  sifter.  The  material  was  sifted  for  7.5  minutes  in  the  sift/pulse  mode  on  a 
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setting  of  8.  Each  sieve  was  removed,  weighed,  tared,  cleaned,  and  weighed  again  to  get 
the  mass  of  material  of  that  fraction.  These  masses  were  then  recorded  and  converted  to 
percent  of  total  sand  distribution: 

%  Si  =  Wi  *  100  (Equation  5) 

Ws 

Where  %  Sj  is  percent  of  sand  fraction  of  a  given  sand  grade  (after  removal  of  organics 
and  fines)  and  Wi  is  the  weight  of  sand  on  a  given  sieve. 

Suspended  Sediment  Concentration 

Suspended  sediment  concentrations  were  measured,  using  the  same  lab 
procedures  as  previous  studies  to  allow  comparisons  and  extension  of  the  sediment 
concentration  data  set  to  the  historical  record  for  1997-1998  at  WI  (Dean  et  al,  1998). 
Samples  were  refrigerated  at  4  °C  until  suspended  sediment  concentrations  were 
measured.  Laboratory  procedures  are  outlined  in  Standard  Methods  for  the  Examination 
of  Water  and  Wastewater  (Eaton  et  al.,  1995).  This  procedure  begins  by  weighing  a  clean 
oven  dried  (at  103  °C)  0.70  um  fiberglass  filter  which  is  placed  on  a  vacuum  filtration 
device.  A  known  volume  of  approximately  75-ml  of  well-stirred  sample  was  extracted 
while  the  sediment  was  in  motion  and  uniformly  distributed.  The  subsample  was  filtered. 
Using  a  vacuum  filtration  device  that  was  thoroughly  rinsed  with  distilled  water. 

Filtration  was  continued  for  an  additional  three  minutes  following  rinsing  to  remove 
excess  water.  The  filter  with  residue  was  removed  and  placed  in  a  numbered  aluminum 
dish,  dried  at  103  °C  for  at  least  three  hours,  removed  to  a  dessicator,  allowed  to  cool,  and 
then  weighed  on  an  analytical  balance.  The  dry  filter  weight  was  subtracted  from  the 
total  dried  residue  and  filter  weight  to  obtain  the  suspended  sediment  weight,  which  was 
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divided  by  the  extracted  sample  volume  to  get  a  suspended  sediment  concentration  in 
mg/1.  Quality  assurances  and  checks  (QA  /  QC)  were  made  with  every  10  to  12  samples 
to  ensure  lab  procedures  are  producing  accurate  measurements  (Eaton,  et  ah,  1995,  Dean 
et  al.,  1998). 

Sediment  flux  was  the  product  of  the  suspended  sediment  concentration  (mg/1) 
and  the  discharge  (1/s)  at  the  time  the  sample  was  taken.  This  product  is  a  sediment  flux 
of  mass  per  unit  time  (mg/s).  This  sediment  flux  was  convert  into  units  of  kg/d.  To 
compare  flux  rates  between  the  W1E  and  WIN  watersheds  the  flux  was  divided  by 
drainage  area  to  yield  unit  sediment  flux  (kg/d/km  ) 

STATISTICAL  METHODS 

Field  samples  represent  a  relatively  small  sample  of  conditions  observed  that  vary 
greatly  in  time  and  space,  so  statistical  relations  are  emphasized  over  individual 
observations.  Extrapolation  of  data  outside  the  observed  ranges  may  result  in  erroneous 
conclusions.  Soil  erosion  models  were  employed  to  assist  with  the  interpretation  of 
observations.  Given  difficulties  with  model  calibration  and  validation,  limited  time,  and 
resources  available  to  this  study,  modeling  was  only  exploratory  and  results  should  be 
considered  preliminary.  Field  measurements  of  erosion,  sediment  concentration,  and 
storage  were  piecemeal  so  some  was  qualitative  by  necessity.  Statistical  test  were  done 
using  data  analysis  tools  for  Microsoft  Excel  for  Windows  95  (Microsoft,  1995)  and  SAS 
procedures  (SAS  Institute,  1998). 

Reeression  Analysis 

Regression  analysis  was  used  to  functionalized  stage  discharge  relationships, 
deposition  along  buffer  zones,  and  analysis  of  suspended  sediment  concentration  and 
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discharge  in  determining  sediment  rating  curves  (Meade,  Yuzyk,  and  Day,  1990;  James 
1998).  Data  were  stratified  by  season,  type  of  precipitation  event,  rising  versus  falling 
limb  of  hydrographs,  year,  and  land-use  change  to  assess  conditions  and  determine 
sediment  load  variability  as  a  function  of  discharges  (Meade,  1982). 

Coefficient  of  Variation  (R2) 

The  coefficient  of  variation  was  used  to  measure  the  percent  variability  of  the 
observed  dependent  variables  as  explained  by  the  variability  of  observed  independent 
variables.  This  measure  was  used  in  regressions  of  flow  stage  on  discharge,  depositional 
mat  unit  weight  distributions  on  distance,  and  sediment  concentrations  (rating  curves)  on 
discharge.  A  high  value  of  R2  indicates  a  strong  association  between  independent  and 
dependent  variables  and  represents  the  percent  variance  in  the  dependent  variable  as 
explained  by  the  independent  variables  (Moore,  1993). 

T-test 

The  t-test  method  is  used  to  determine  if  the  difference  between  two  means  is 
significant.  This  test  was  used  on  Hypothesis  one,  two,  three,  and  on  quality  assurances 
and  checks  on  laboratory  procedures.  For  each  Hypothesis  the  null  hypothesis  (Ho)  was 
evaluated;  that  is  that  there  was  no  difference  between  means.  Ho  was  rejected  if  the  t-test 
indicated  significant  difference  at  the  alpha  =  0.05  level.  Null  Hypothesis  1  was  tested  by 
the  difference  between  observed  soil  loss  from  dirt  roads  against  forested  control  sites. 
Null  Hypothesis  2  was  tested  by  the  difference  between  sediment  deposition  and  particle 
size  distributions  at  various  positions  in  the  buffer  zones.  Null  Hypothesis  3  was  tested 
by  the  difference  between  storm  unit  sediment  flux  (kg/d/km2)  at  W1E  and  WIN  for  five 
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storm  events.  Quality  assurances  and  checks  were  paired  repetitive  samples,  which  were 
tested  for  significant  difference  from  zero. 

ANCOVA 

Analysis  of  covariance  (ANCOVA)  evaluates  quantitative  variables  by 
classification  with  a  qualitative  variable  by  calculating  an  F-statistic  which,  when 
compared  to  an  F-value  based  on  the  alpha  level  and  degrees  of  freedom,  test  the 
significance  of  differences  in  regression  slopes  under  the  stratification  of  the  descriptive 
qualitative  variables.  Using  analysis  of  covariance  (ANCOVA)  inferences  on  the 
differences  between  previous  study  data  and  this  study  data  to  determine  compatibility. 
ANCOVA  was  used  to  evaluated  WIN  and  W1E  unit  flux  unit  flow  rating  curves  and  a 
pre-timber  harvesting  curve  against  a  post-timber  harvesting  event  to  reflect  any 
significant  difference. 


CHAPTER  III 


DATA  ANALYSIS  AND  RESULTS 

SAMPLED  EVENTS 

Nine  storm  events  were  sampled  during  from  February  1999  to  January  2000. 

This  represented  most  of  the  larger  storms  during  the  study  period  plus  a  few  small  events 
(Figure  3-1).  Anticipated  limitations  of  this  study  included  the  need  to  reach  the  study 
area  quickly  before  the  onset  of  storms  to  capture  the  initial  rising  limb  response.  Travel 
time,  precipitation  spatial  variability,  and  distances  between  sample  sites  constrained  the 
collection  of  data.  During  storms,  priority  was  given  to  suspended  sediment  sampling  but 
field  observations  across  the  basin  were  also  necessary  to  evaluate  the  spatial  nature  of 
responses  from  sediment  sources,  runoff,  buffer  zone  effectiveness,  and  wetland 
inundation.  Total  rainfall  for  the  study  period  (February  1999  to  December  1999)  was 
61.5cm  compared  with  a  total  of  74  cm  which  was  about  40  cm  below  the  annual  mean 
(unpublished  data,  Carbone,  2000;  US  Dept,  of  Commerce,  1992). 
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9B  Precip  1999 


Date-Time 


Figure  3-1.  Sampled  Events  and  the  precipitation  record  at  Bravo  9  Rain  Gage 
(Rain  Data  Unpublished,  Carbone,  2000) 

SEDIMENT  EROSION 
Rill  Erosion 

Rill  erosion  susceptibility  along  the  dirt  roads  was  evaluated  by  monitoring  cross- 
section  surveys  along  the  hillslope  at  R3  and  a  gully  at  site  R2  to  estimate  profile  changes 
through  time  and  cross-sections  and  volumes  were  calculated  from  net  changes  between 
surveys.  The  R3  hillslope  had  a  profile  length  of  98  m,  width  of  4  m,  and  an  average 
slope  of  5%.  The  R2  gully  survey  site  was  2.6  m  long,  2.3  m  wide,  with  a  slope  of  1%. 
R3  Hillslope  Profile 

The  R3  profile  was  measured  in  two  ways.  Four  separate  cross-sections  and  the 
right  center  (1.5  meters  from  the  right  side,  viewing  uphill)  profile  were  monitored  from 
June  1999  to  November  1999.  This  represented  the  areas  of  maximum  erosion  activity. 
During  the  six  month  period  (mid  June  1999  to  early  November  1999)  the  cross-section 
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soil  loss  was  15.68  kg/m2.  This  represents  erosion  of  approximately  140  tons/ac/yr 
(Table  3-1)  based  on  a  road  surface  area  calculated  as  the  segment  length  (m)  times  the 
roadway  mean  width  (4  m).  If  this  rate  continued  for  another  seven  months  the  soil  loss 
would  exceed  the  tolerance  (T)  factor  in  the  county  soil  survey  of  3  tons/ac/yr 
(Lawerance,  1978).  The  greatest  soil  loss  was  measured  near  the  base  of  the  hillslope. 
This  indicates  that  sediment  is  being  effectively  transported  downslope  and  off  the  slope. 
Local  deposition  began  to  occur  beyond  the  98  m  point  of  the  hillslope  base  and  extended 
about  five  meters  until  it  reached  the  forest  edge.  At  this  location  the  road  gradient  is 
minimal  at  one  to  two  percent.  The  gradient  slopes  away  from  the  hillslope  toward  the 
WIN  tributary.  Field  observations  monitored  the  continued  existence  of  a  sediment 
plume  that  is  formed  at  the  edge  of  the  forest  and  extended  eight  meters  below  the  road 
(Figure  3-6).  Prior  to  the  study,  most  dirt  roads  were  graded  leaving  the  roadway  at  R3 
smooth  with  loose  sediment  material.  Subsequent  surveys  characterize  the  development 
of  rills  mainly  on  the  lower  south  side  of  the  road  looking  uphill  (Figures  3-2  to  3-5;  cf. 
Figure  2-7).  Each  cross-section  had  periods  of  deposition  and  erosion  reflecting  the 
downward  movement  of  sediment.  The  dominant  areas  of  erosion  and  rill  development 
were  on  the  lowest  side  of  the  road  where  flow  velocities  are  maximized.  Another  area  of 
high  rill  development  was  the  embankment  of  spoils  from  grading  on  both  sides  of  the 
roadway.  Deposition  occurred  mainly  on  the  north  side  at  the  base  of  the  embankments 
where  the  roadway  had  the  smallest  gradient. 
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R3  Cross-sections 

0m-15m  15m 

-  40m  40m 

-75m 

Average  Area  of  Net  Erosion(mA2) 

0.0352 

0.0501 

0.0655 

Segment  Length  (m) 

15 

25 

35 

Volume  (L) 

528.0 

1252.8 

2290.8 

Bulk  Density  (g/cc) 

1.52 

1.55 

1.59 

Net  erosion  (kg) 

803 

1942 

3642 

Unit  Net  Erosion  (kg/mA2) 

10.70 

15.53 

20.81 

Average  Unit  Erosion  (kg/m2) 

15.68 

Table  3-1.  R3  Cross-section  Soil  Losses  (Distance  Measured  from  top  of  Hillslope) 

R3  Cross-section  (0  m) 

N  Net  Area  Lost  (-.003  m2)  S 


4  Jun  99 
19-Nov-99 


Width  (m) 


Figure  3-2.  R3  Cross-section  Soil  Losses  at  0  m  (See  page  22) 


Figure  3-3.  R3  Cross-section  Soil  Losses  at  15  m 


Figure  3-4.  R3  Cross-section  Soil  Losses  at  40  m 
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Figure  3-6.  Sediment  Deposition  below  R3  Hillslope  Base,  Looking  Downhill 
Toward  WIN  Channel  Below  Earthen  Dam 
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The  right  centerline  profile  indicated  a  continual  loss  of  sediment  down-slope 
(Figure  3-8).  This  profile  also  indicates  efficient  sediment  transport  with  no  deposition. 


Figure  3-7.  R3  Right  Center-line  Profile  Survey 

i 

The  WEPP  Hillslope  model  was  run  with  a  one-year  simulation  to  simulate  soil 
erosion  and  to  compare  results  with  the  field  measurements  of  the  R3  profile.  These 
preliminary  modeling  results  are  intended  to  analyze  a  general  trend  for  this  hillslope 
using  basic  model  parameters  of  management  practices,  slope  profile,  soil  profile,  and 
regional  climate  generators.  The  model  results  were  similar  to  the  observed  general  trend 
of  detachment  and  profile  lowering  along  the  entire  length  with  no  deposition  (Figure  3.9) 
but  modeled  a  lower  rate  of  erosion  at  0.85  kg/m  /yr  or  0.35  kg/m  /(5  months)  as 
compared  to  the  R3  cross-section  soil  loss  rate  (15.68  kg/m2/(5  months)  (Table  3-2). 
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R3  WEPP  Results 

Management 

Range-Graze 

Soil  Name 

Vaucluse 

Average  Annual  Precip(mm) 

1160.90 

Average  Annual  Runoff  (nun) 

260.30 

Average  Annual  Soil  Loss  (kg/m 2) 

0.85 

Detachment  Length  (m) 

98.00 

Average  Deposition  (t/ha) 

0.00 

Deposition  Length  (m) 

0.00 

Table  3-2.  R3  Profile  WEPP  Model  Results 


Limitations  of  this  simulation  include  the  lack  of  available  management 
parameters  for  the  WEPP  model  for  Mill  Creek  conditions.  For  example,  there  are  no 
parameters  for  road  surfaces  or  the  Vaucluse  soil  profile.  Use  of  this  model  in  forested 
areas  with  timbering  practices  therefore,  required  some  generalization  and  manipulation 
of  available  parameters  in  this  area.  Past  WEPP  model  validations  have  indicated  that 
rangeland  interrill  erodibility  values  were  underpredicted  and  rill  erodibility  values  were 
overpredicted.  Although  these  errors  were  not  as  much  as  with  cropland  values  used  on 
recently  graded  forest  roads,  caution  is  advised  when  utilizing  this  model  outside 
agricultural  areas  (Elliot  et  al.,  1995).  A  bare  rangeland  management  practice  was 
constructed  to  match  the  dirt  road  profile.  A  Vaucluse  soil  profile  was  built  from  field 
measurements  to  reflect  the  change  caused  by  grading  in  the  county  Vaucluse  soil  profile. 
The  field-surveyed  length-slope  values  of  the  profile  and  the  WEPP  regional  climate  data 
from  Aiken,  SC  were  used  in  the  simulation.  This  model  will  require  more  in-depth 
parameterization  and  field  analysis  to  better  represent  the  local  physical  properties  and 
management  practices  if  used  as  a  prediction  and  management  tool.  Given  the 
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substitution  of  surrogate  climate,  land-use,  and  soil  parameters,  greater  credence  is  given 


Figure  3-8.  Preliminary  WEPP  Hillslope  Profile  Model  Graphical  Results.  Top 
curves:  Initial  Profile  Being  the  Solid  Line  and  the  Predicted  Profile  as  the  Dotted 
Line  and  Bottom  Profile  predicting  the  Soil  Loss  Rate  along  the  Profile 


The  WEPP  hillslope  model  outputs  a  display  of  the  dirt  road  initial  profile  with 
the  solid  line  and  the  resulting  profile  after  the  one-year  simulation  with  the  dotted  line 
(Figure  3-9).  The  trend  of  continuous  soil  loss  is  reflected  in  the  model  with  the 
reduction  in  elevation  along  the  profile.  The  bottom  portion  of  the  soil  loss  graph 
indicates  the  change  in  erosional  rates  along  the  profile.  The  maximum  rate  occurred 
around  95  m  down-slope.  The  rapid  decline  of  the  erosion  rate  from  the  maximum  point 
to  the  segment  end  at  98  m  implies  that  deposition  may  start  to  occur  soon  after  this  point. 


43 


This  would  be  consistent  with  field  observations  at  the  site.  Beyond  the  98m  point  of  this 
profile,  the  hillslope  changes  gradient  as  it  intersects  two  other  dirt  roads  and  is  only  three 
meters  from  the  forest  woodline.  The  WEPP  corroborates  the  lower  erosion  rates 
observed  at  the  uppermost  cross-section  where  there  is  less  runoff.  The  highest  erosion 
rates  occur  near  the  slope  base  at  the  greater  slope  lengths  where  the  runoff  reaches 
maximum  volumes  and  velocities  and  has  the  greatest  detachment  potential.  Beyond  96 
m  the  breaks  in  slope  and  increased  infiltration  into  recent  sandy  deposits  decrease  the 
sediment  transport  capacity  and  result  in  further  deposition. 

R2  Gully  Profile 

The  R2  gully  erosion  site  was  different  than  the  R3  hillslope  site  in  that  was  an 
established  rill  on  which  road  runoff  was  actively  maintained  (Figure  2-8).  There  was  no 
dominant  erosion-susceptible  area  within  the  gully  and  for  the  most  part  observed 
sediment  movement  across  the  channel  was  uniform.  The  R2  segment  five-month  soil 
loss  was  1.42  kg/m2  based  on  the  product  of  the  soil  loss  area  (0.21  m2)  and  the  distance 
between  cross-sections  (2.6  m)  divided  by  the  segment  area  5.89  m2.  This  soil  loss  is 
greater  than  that  of  R3  and  is  representative  of  the  erosive  potential  of  rill  erosion  over 
loose  sediment.  As  with  the  R3  hillslope,  the  R2  gully  had  a  change  of  gradient  as  it  left 
the  road  and  entered  the  forest  wood-line.  Unlike  the  R3  hillslope  site,  the  R2  gradient 
increased  to  6  percent.  This  resulted  in  a  long  sediment  plume  that  had  a  length  of  78.5 
meters  as  measured  from  the  forest  edge. 

Interill  Small  Box  Plots 

The  small  box  plot  sediment  from  the  2-foot  by  2-foot  wooden  boxes  indicate  soil 
susceptibility  to  entrainment  forces  by  raindrop  splash  and  sheet  erosion.  Evidence  of 
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splash  erosion  could  be  seen  on  the  sides  of  the  box  with  negligible  sediment  on  the  tops 
of  board  walls.  This  indicates  that  little  sediment  is  being  splashed  into  or  out  of  the  box 
but  could  represent  some  losses  or  at  least  a  lag  in  sediment  reaching  the  box  drain  spout. 
This  error,  described  in  Lai  (1994)  as  the  edge  effect,  was  minimized  by  ensuring  similar 
ground  conditions  surrounded  the  boxes.  Eleven  box  sample  plots  were  maintained  from 
mid-June  1999  to  mid-September  1999  for  five  storm  events.  Increased  vehicular  traffic 
and  human  activity  in  the  study  area  compromised  the  boxes  after  September  1999.  Five 
control  boxes  were  established  in  representative  forested  areas  and  six  test  boxes  were 
established  on  dirt  roads.  The  dirt  road  unit  sediment  loss  (kg/m2)  means  were 
hypothesized  to  be  significantly  higher  than  the  forested  control  group.  The  dirt  road 
sediment  means  were  tested  and  found  to  be  significantly  greater  than  forested  sediment 
at  the  0.05  alpha  level  (t-Stat  =  2.3152  and  P  (T<t)  =  0.0342)  (Figure  3-2). 


Interill  Box  Plots 


Forested  Dirt  Road 

Box  Type 


Figure  3-9.  Interrill  Box  Plot  Mean  Unit  Mass 
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DEPOSITION 


Sediment  deposition  was  noted  at  several  locations  within  the  study  area. 
Deposition  was  observed  mainly  at  locations  where  road  runoff  entered  the  forest  edge. 
Additionally,  deposition  occurred  on  road  slopes  where  decreases  in  gradient  reduce  flow 
energies.  Two  other  areas  of  deposition  were  within  the  wetland  and  along  channel  beds. 
In  the  wetlands,  evidence  of  fines  coating  the  litter  layer  was  observed  as  a  result  of 
inundation  and  overbank  flow.  Within  the  ditched  channels  of  W1E  and  WIN,  some 
temporary  sediment  deposition  occurred  but  for  the  most  part  the  channels  remained  clear 
as  the  thin  veneer  of  bed  sediment  was  removed  with  subsequent  storms.  One  exception 
was  the  W1E  tributary  draining  the  BRM-19  rifle  range,  which  remained  filled  with 
sandy  sediment.  Field  observations  of  changes  in  morphology  of  the  channel  sand 
deposits  indicate  these  bed  materials  were  moving  during  storms.  Another  exception  is 
on  the  W1E  channel  below  the  confluence  with  this  BRM-19  tributary,  where  the  ditch 
morphology  is  observed  by  channel  aggradation.  The  wide  shallow  channel  in  this  reach 
indicates  sediment  storage  in  the  main  channel. 

Considerable  sediment  storage  also  occurs  below  roads  in  the  transition  zone  at 
the  wetland  margins.  Buffer  zone  deposition  efficiency  was  evaluated  by  measuring  the 
reduction  of  sediment  being  transported  through  the  transition  zone  between  sediment 
sources  (a  gully)  at  the  road  edge  and  the  wetlands  toward  the  channel.  This  transition 
zone  acts  as  a  buffer  zone  as  flow  velocities  are  reduced  upon  entering  the  gentle  slopes 
of  the  vegetated  wetland.  The  unit  weight  (kg/m2)  of  sediment  deposition  was  measured 
on  five  occasions  over  a  six-month  period.  Transects  were  sampled  at  five  different 
locations  with  less  than  a  3  percent  slope  (Figure  2-9).  Each,  sampling  consisted  of  four 
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mats  along  a  transect  at  ten  meter  intervals.  The  unit  weights  along  each  of  the  five  mat 
transects  indicate  that  a  large  portion  of  the  sediment  was  trapped  in  the  first  10  m  of  the 
buffer  zone  with  a  exponential  decline  reaching  a  minimum  around  30  m  and  an  slight 
increase  in  deposition  at  the  40  m  width  (Figure  3-10).  This  indicates  that  the  majority  of 
deposition  occurred  in  the  first  portion  of  the  buffer  zone  but  occasional  flooding  near  the 
channel  presumably  caused  the  slight  increase  at  the  40  m  position  by  overbank 
suspended  sediment  or  slackwater  conditions.  This  sedimentation  from  slackwater  or 
overbank  flow  may  not  be  applicable  on  hillslope  transect  away  from  channels  or 
wetlands.  The  mean  unit  mass  deposited  at  the  first  10  m  buffer  width  was  significantly 
greater  than  the  combined  mean  deposition  weights  at  the  other  three  sites  at  the  0.05 
alpha  level  (t-Stat  =  2.27  and  P  (T<=t)  is  0.04).  Thus,  buffer  widths  from  20  m  to  40  m 
would  be  far  less  effective  at  trapping  sediment  than  the  narrow  strips,  based  on  these 
samples.  Increased  runoff  during  larger  flows,  less  vegetation,  and  areas  of  steeper  slopes 
may  have  different  results. 
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Deposition  Mats 


Figure  3-10.  Sediment  Deposition  on  Mats  Occurred  mostly  Near  the 
Edge  of  the  Wetland.  Each  Storm  was  Measured  at  a  Different  Transect 

of  Four  Mats 


A  statistical  regression  of  the  unit  mass  against  transect  distance  explains  70%  of 
the  variance  in  sediment  deposition  (Figure  3-11).  At  the  30  m  buffer  width,  the 
statistical  model  predicts  sediment  deposition  is  essentially  zero  (-0.02).  Due  to  increased 
sedimentation  caused  by  inundation  near  main  channels  in  the  wetlands,  the  40  m  buffer 
width  may  not  be  necessary  in  this  environment  for  this  size  runoff  events  if  an 
insignificant  amount  of  sediment  is  being  transported  past  the  30  m  width.  However,  if 
continued  long-term  deposition  causes  steeper  gradients  that  may  result  in  sediment  being 
transported  further  into  the  buffer  zone. 
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Deposition  Mats 


Figure  3-11.  Sediment  Deposition  as  a  Function  of  Distance  along  a  Transect 


Deposition  mat  particle  textural  analysis  reinforced  the  unit  mass  findings. 
Statistical  averages  of  loss-on-ignition  (LOI)  and  sediment  textures  from  wet  sieving  for 
twenty  mat  samples  are  shown  in  Figure  3-12  as  percent  weight.  These  sedimentologic 
results  show  that  the  greatest  percentage  of  sand  is  deposited  in  the  10  m  buffer  width 
along  with  the  least  percentage  of  fines  or  organics.  This  indicates  that  the  larger  sand 
particles  are  settling  out  quickly  with  the  increased  roughness  and  decreased  slope.  The 
fine  sediment  is  being  transported  further  along  the  transect.  There  is  also  an  increase  in 
organic  matter  as  the  runoff  mobilizes  the  litter  layer.  The  little  change  in  texture  or 
organics  between  the  30  m  and  40  m  buffer  widths  may  reflect  deposition  dominated  by 
over-bank  inundation  near  the  main  channel  of  the  east  tributary. 


49 


Deposition  Mat  Sediment  Textures  and  Organics 
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Figure  3-12.  Deposition  Mat  LOI  and  Wet-Sieving  Results  (Percentages  by  Weight) 
Values  at  each  Horizontal  Position  are  an  Average  of  Five  Events  at  that  Site 

Sand  distributions  also  reflect  deposition  along  the  transect.  Sand  textures 
measured  by  sonic-sieving  are  shown  in  Figure  3-13.  The  coarser  sands  (>  0.5mm) 
generally  decrease  with  an  increase  in  buffer  width.  This  fining  away  from  the  source  is 
also  reflected  in  the  distribution  of  the  fine  and  very  fine  sand  particles  (<0.25mm  and 
>0.63mm).  There  is  a  sharp  increase  in  the  fine  and  very  fine  sand  particles  with  the 
increase  in  buffer  width  reaching  a  peak  at  the  30  m  width.  Typically,  this  distribution 
indicates  that  deposition  is  occurring  in  the  lower  velocity  flow  environments  after  the 
larger  particles  had  settled.  Very  little  material  coarser  than  sand  (>  2mm)  was  being 
transported.  One  large  pebble  was  found  in  a  30  m  sample  and  may  indicate  scour  of  an 
erosional  lag  from  an  adjacent  surface.  Presence  of  coarser  particles  may  reflect 
substantial  distances  of  sediment  transport  in  large  storm  events.  Therefore,  these  results 
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conclude  that  sediment  deposition  is  substantially  reduced  with  distance  traveled  across  a 


buffer  zone  transect. 


Figure  3-13.  Sand  Textures  of  Sediment  on  Deposition  Mat  after  Removal  of 

Fines  and  Organics 

SEDIMENT  YIELD 

This  section  addresses  several  factors  that  contribute  to  sediment  yield.  Topics 
include  discharge,  suspended  sediment  concentrations,  suspended  sediment  storm  flux, 
hysteresis,  and  study  area  land-use  changes. 

Discharge 

Nine  discharge  readings  were  taken  at  gage  7  on  W1E  in  order  to  establish  a 
stage-discharge  relationship  and  functionalize  discharge  (Figure  2-3  and  Table  2-2).  A 
second-order  polynomial  had  a  high  explained  variance  (R2  =  0.98).  Residual  analysis 
had  no  indication  of  change  through  time  for  the  discharge  measurement  period  of  June 
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through  September.  Statistical  inferences  utilizing  a  confidence  interval  of  95%  indicates 
an  accuracy  of  +  1.59  (1/s).  Manual  staff  gage  readings  and  continuous  stage  recorders 
monitored  discharge.  A  statistical  regression  of  manual  staff  gage  readings  functionalized 
the  continuous  fifteen-minute  recordings  of  mV  values  provided  an  expression  for  stage 
as  a  function  of  logger  output  (R2of  0.95). 

The  variance  was  a  result  of  limited  staff  gage  readings,  the  +  1  cm  margin  of 
error  in  the  pressure  transducers  used,  and  the  variance  caused  by  the  fifteen  minute  time 
interval.  The  function  allowed  the  conversion  of  logger  mV  values  to  gage  heights  values 
and  subsequently  into  discharge  values  by  applying  the  stage  discharge  function  in 
(Figure  2-3).  Equipment  failures  resulted  in  a  period  of  lost  data  prior  to  June  and  after 
late  December  1999. 

A  preliminary  continuous  discharge  record  for  this  period  was  obtained  and 
utilized  in  examining  the  discharge  data  prior  to  and  after  the  visual  staff  readings  to  gain 
general  insight  into  storm  events  response.  Since  further  data  collection  is  needed  to  gain 
a  better  relationship  between  continuous  logger  readings,  staff  gage  heights,  and 
discharge,  greater  credence  was  given  to  manual  staff  gage  readings  in  the  calculation  of 
discharge.  Instantaneous  discharge  calculations  for  corresponding  suspended  sediment 
samples  allowed  further  application  with  suspended  sediment  flux.  Hydrographs  were 
constructed  to  depict  the  arrival  and  timing  of  discharge  peaks  and  allow  calculation  of 
storm  discharge  (e.g.  Figure  3-14). 
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Figure  3-14.  Single  Storm  Runoff  Hydrograph  for  12  February  1999  with 
Suspended  Sediment  Flux  at  the  three  Gages.  Note  the  Early  Arrival  of  high 

Sediment  Fluxes  at  WIN  and  W1 

Base  flow  was  separated  from  storm  flow  for  five  storm  events  by  applying  a 
method  for  basins  less  than  20  square  miles.  A  straight  line  was  drawn  from  the  initial 
point  of  rise  at  a  rate  of  0.05  cfs  per  square  mile  per  hour  (Hewlett  and  Hibbert,  1967). 
Instantaneous  storm  flows  were  calculated  by  subtracting  base  flow  from  the  total 
discharge.  For  the  total  storm  runoff,  the  area  under  the  curves  for  Wl,  W1E,  and  WIN 
was  calculated  as  the  sum  of  the  products  of  instantaneous  flows  and  their  respective 
duration  (Figure  3-14). 
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Suspended  Sediment  Flux 


The  calculated  total  suspended  sediment  flux  for  an  individual  storm  provides  the 
sediment  yield  for  that  event  at  the  three  gage  sites  to  test  for  differences  between  W1E 
and  WIN.  Total  sediment  storm  flux  and  discharge  magnitudes  were  calculated  for  five 
1999  events  (12  Feb  99, 16  Jun  99,  19  Jul  99,  4  Oct  99,  and  10  Jan  00)  (Table  3-3).  Due 
to  differences  between  1997  and  1998  sediment  concentrations  discussed  later,  no  attempt 
was  made  to  compare  sediment  fluxes  between  WIN  and  W1E  for  the  earlier  period. 


Total  Q  (m3) 

Total  Storm  Sediment  Flux  (kg) 
W1  W1E  WIN 

Date 

11.4 

5.2 

3.53 

1.68 

12-Feb-99 

78.5 

117.2 

47.7 

69.5 

16-Jun-99 

4.3 

19.2 

1 

18.2 

19-Jul-99 

9.9 

0.068 

0.679 

2.19 

4-Oct-99 

3207.8 

633.2 

392.7 

240.5 

10-Jan-00 

Table  3-3.  Selected  Storm  Flux  and  Discharge  Totals 


Field  observations  indicated  sediment  concentrations  came  in  pulses  at  varying 
time  intervals  from  the  two  tributaries.  WIN  tended  to  peak  ten  to  fifteen  minutes  from 
the  centroid  of  precipitation.  In  response  to  quick  intense  storms  with  low  discharge,  like 
12  July  1999,  WIN  peaked  at  five  minutes  with  a  peak  suspended  sediment  concentration 
of  2000  mg/1.  W1E  had  little  sediment  response  to  this  storm  with  a  peak  suspended 
sediment  concentration  of  only  71  mg/1.  The  maximum  rainfall  intensity  for  19  July  1999 
was  0.84  mm/minute  and  a  total  precipitation  of  1 1.2  mm.  With  a  low-intensity  storm 
such  as  4  October  1999,  WIN  responded  with  a  high  suspended  sediment  peak 
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concentration  of  371  mg/1  about  fifteen  minutes  after  the  precipitation  centroid.  During 
this  low-intensity  storm  event,  the  sediment  pulse  in  W1E  was  delayed  about  45  minutes 
past  the  precipitation  centroid  and  had  a  maximum  suspended  sediment  concentration  of 
21  mg/1.  The  main  difference  in  the  storm  events  was  that  the  later  storm  had  a  lower 
rainfall  intensity  (0.25  mm  per  minute)  than  the  first  storm.  The  sediment  concentration 
differences  between  these  two  storms  reflect  the  erosive  potential  of  high  intensity  storm 
events. 

Sediment  contributions  from  the  W IE  or  WIN  tributaries  do  not  provide  constant 
proportions  of  the  total  sediment  load.  Generally,  WIN  is  the  main  contributor  to  storm 
flux  in  short  intense  convectional  storms.  In  the  longer  duration  frontal  storms,  W1E 
contributed  the  most  sediment.  The  mean  sediment  fluxes  for  these  five  storm  events 
were  tested  to  see  if  there  was  significant  difference  between  the  two  tributaries.  No 
significant  difference  in  sediment  flux  was  found  between  the  two  channels  at  the  0.05 
alpha  with  (t-Stat  =  0.69  and  P  (T<=t)  =  0.53).  Since  WIN  tends  to  contribute  high 
amounts  of  sediment  during  the  shorter  events,  the  small  storm  means  were  tested 
separately  and  again  there  was  no  significant  difference  at  the  0.05  alpha  level  with  (t-Stat 
=  1.67  and  P(T<=t)  =  0.19). 

The  assumption  of  equation  two,  (Qswi  =  Qswin  +  Qswie),  is  that  the  sediment 
discharge  of  the  main  channel  is  the  sum  of  the  suspended  sediment  discharge  of  the  two 
tributaryes  immediately  upstream.  This  assumption  was  tested  using  a  t-test  on  the 
sediment  discharges  for  W 1  verses  WIN  and  W1E.  Using  samples  taken  at  each  of  the 
three  sites  at  no  more  than  three  minutes  apart,  the  sediment  flux  of  the  main  channel  was 
found  to  have  no  significant  difference  from  the  sum  of  the  north  and  east  channels  at  the 
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0.05  alpha  level  (t-stat  =  0.80;  P  =  0.43).  This  test  confirms  that  contributing  discharge 
and  suspended  sediment  load  at  WIN  can  be  calculated  with  reasonable  accuracy  by 
measuring  water  and  sediment  discharges  in  the  main  channel  and  the  east  tributary. 

To  chart  a  general  relationship  between  sediment  flux  and  discharge  the  storm 
flux  data  were  plotted  against  the  log  of  storm  discharge  (Figure  3-15).  The  results 
indicate  a  general  rise  of  storm  flux  with  storm  discharge  but  considerable  scatter.  Storm 
flux  and  storm  discharge  was  also  evaluated  against  maximum  precipitation  intensity  and 
total  rainfall  depth.  Sediment  load  varied  erratically  with  rainfall  intensity  but  there  is  the 
suggestion  that  rainfall  intensities  are  associated  with  higher  sediment  fluxes  in  the  north 
tributary  than  in  the  east  tributary  (Figure  3-16).  A  sharp  increase  in  sediment  flux  occurs 
with  total  rainfall  suggesting,  as  expected,  that  more  sediment  is  moved  by  larger  events 
(Figure  3-17).  The  high  variability  in  sediment  flux  with  discharge,  precipitation 
intensity,  and  precipitation  totals  call  for  a  more  thorough  examination  of  the  relationship 
between  sediment  concentration  and  runoff. 
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Total  Stom  Discharge  (mA3) 


Figure  3-15.  Storm  sediment  Flux  and  Storm  Discharge  Totals  (Same  storm  events 

in  Table  3-3) 


Figure  3-16.  Storm  Flux  Response  to  Max  Rainfall  Intensity  (Same  storm  events  in 


Table  3-3) 


Small  Event  Sediment  Flux  Response 


Figure  3-17.  Storm  Flux  Responses  to  total  Precipitation  (Same  storm  events  in 

Table  3-3) 


Suspended  Sediment  Concentrations  and  Rating  Curves 

Calculated  as  the  product  of  measured  suspended  sediment  concentrations  of  Wl, 
WIN,  and  W1E  and  their  respective  flow  discharges.  Suspended  sediment 
concentrations  were  highly  variable  due  to  fluctuations  in  the  timing  of  pulses  of 
discharge  and  sediment  coming  from  each  tributary.  Suspended  sediment  concentration 
data  tables  are  found  in  Tables  D1  through  D12,  (Appendix  D).  Basic  suspended 
sediment  concentration  statistics  reflect  the  wide  range  of  concentrations  (Table  3-4). 
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Gage  Site 

Min 

Max 

Mean 

Std  Dev 

□O 

WIN  (G6) 

6.10 

2529.41 

284.03 

W1E(G7) 

6.33 

394.74 

83.51 

W1  (G1) 

6.54 

1181.82 

126.67 

172.78 

60  | 

Table  3-4.  Basic  Statistics  of  1999  Suspended  Sediment  Concentration  (mg/1) 


by  Sampling  Site 


In  order  to  expand  this  study’s  results  with  other  studies  of  the  watershed,  several 


suspended  sediment-rating  curves  were  developed  to  analyze  sediment  response  to 
discharge  under  different  conditions.  These  rating  curves  relate  instantaneous  (single 
sample)  sediment  concentration  measurements  to  the  corresponding  streamflow  delivered 
from  the  stage  observations.  Both  suspended  sediment  concentrations  and  unit  flux  rating 
curves  allow  comparison  between  the  two  subbasins  of  the  north  and  east  tributaries. 
Stratification  of  data  was  necessary  to  gain  insight  into  sediment  responses.  Extensive 
land-use  changes  and  an  unusually  low  rainfall  allowed  examination  of  responses  to  these 


changes. 

Comparison  with  Previous  Studies 

A  large  number  of  suspended  sediment  samples  were  collected  and  processed 
between  1997  and  1998  (Dean  et  al,  1998;  2000).  Analysis  of  covariance  with  the  W1 
data  between  1997-98  and  1999  indicate  significant  difference  is  evident  (Figure  3-18). 
Statistically,  the  1999  data  have  significantly  higher  sediment  concentrations  than  the 
earlier  data  at  the  0.05  alpha  level  with  a  (F=  26.48  and  P  (>F  =  0.0001).  This  may  be  a 
result  of  changes  in  land  use,  particularly  the  increased  logging  activity  and  road  work,  or 
drought  conditions  which  allow  an  increase  in  loose  sediment  particles  to  accumulate 
between  storm  events.  Due  to  this  significant  difference  in  sediment  rating  curves,  the 
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1997  and  1998  data  were  not  included  in  further  analysis  of  the  W1  sediment  data 
collected  by  this  study. 


Changes  in  SS  Discharge  Rating  Curves 


Figure  3-18.  W1  suspended  Sediment  Concentration  Data  from  1997-98-99  as 

Functions  of  Discharge. 


Stratified  1999  Sediment  Rating  Curves 

The  1999  sediment  concentration  data  were  split  into  W1  and  W1E  data  series. 
This  allowed  for  comparisons  of  the  two  tributaries  with  the  inference  that  any 
differences  resulted  from  the  WIN  tributary  contributions  to  Wl.  The  total  W1  sediment 
rating  curves  show  a  high  degree  of  variability  in  sediment  concentration  with  discharge 
(Figure  3-19).  Although  regression  equations  are  presented,  the  high  unexplained 
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variance  limits  their  utility  for  long-term  sediment  flux  calculations  and  calls  for 
stratification  of  the  1999  data  (Figure  3-19). 

Several  approaches  were  taken  to  systematically  examine  the  data  for  possible 
reasons  for  the  high  variance  and  improve  the  understanding  of  the  W1  watershed 
sediment  responses  to  discharge.  The  first  approach  was  to  separate  data  based  on  rising 
verses  falling  limbs  of  the  hydrograph.  The  peak  discharge  reading  was  included  in  the 
rising  limb.  The  next  stratification  was  a  seasonal  separation  of  the  data.  In  this 
watershed,  the  storms  were  mainly  convectional  events  from  April  through  September 
and  were  more  typically  frontal  storms  from  October  through  May.  There  were  no  major 
precipitation  contributions  from  hurricane  events  during  this  study. 

The  rising  versus  falling  limb  stratification  resulted  in  an  improved  relation  in  the 
falling  limb  rating  curve  (R2  =0.44)  but  had  a  higher  variability  (R2  =  0. 1 1)  in  the  rising 
limb  rating  curves  especially  in  the  W1  data  as  a  result  of  WIN  contributions.  The  arrival 
of  the  first  flush  of  suspended  sediment  in  WIN  explains  the  weak  rising  limb 
relationships. 

The  quick  response  times  to  runoff  in  this  small  watershed  commonly  results  in 
multiple  pulses  of  sediment,  an  early  pulse  (first  flush),  and  bimodal  peaks  due  to 
variability  of  rainfall  intensities  during  storm  events.  All  of  these  factors  are  present  in 
the  study  area.  Response  times  vary  between  the  W1E  and  the  W1E  tributaryes 
complicating  the  suspended  sediment  concentration  in  the  W1  main  channel. 

The  WIN  site  is  greatly  influenced  by  the  first  flush  phenomenon.  The  peak 
sediment  concentrations  arrive  very  early  within  five  to  ten  minutes  of  the  centroid  of 
precipitation  and  with  intense  storms  within  a  five  to  ten  minutes  of  the  start  of  the 
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precipitation.  In  this  initial  period,  the  concentration  is  spiked  with  a  pulse  of  suspended 
sediment  moving  down  the  WIN  tributary.  The  concentration  peaks  early  in  the  rising 
limb  and  then  begins  to  decrease  as  the  discharge  continues  to  rise.  During  longer  events, 
there  may  be  a  second  pulse  from  WIN  as  sediment  from  further  upstream  arrives  and  the 
first  flush  sediment  has  been  removed.  The  second  peak  is  high  as  the  first  flush.  As  a 
result  of  the  first  flush  phenomenon  W1  suspended  sediment  concentrations  are  highly 
variable  with  discharge  (R2  =  0.1 1)  during  the  first  fifteen  minutes  from  the  centriod  of 
the  precipitation. 

Analysis  of  precipitation  data  and  field  observations  suspended  sediment 
concentration  arrivals  reflected  this  quick  response  from  W IN.  The  fast  increase  of 
suspended  sediment  came  from  direct  road  runoff  that  enters  W IN  about  50m  upstream 
from  the  sampling  site  (Dean  et  al.,  1998;  2000).  To  isolate  the  random  error  associated 
with  this  first  flush,  the  first  fifteen  minutes  of  each  storm  (measured  from  the  centroid  of 
precipitation)  was  separated  from  other  sediment  concentration  measurements.  After 
removal  of  these  early  data  points,  the  variance  in  sediment  concentration  as  explained  by 
discharge  increased  substantially  to  R2  =  0.37  (Figure  3-20).  Comparisons  of  figures  31- 
19  and  3-20  reveal  that  removal  of  the  first  fifeteen  minutes  of  the  hydrographs  had  little 
effect  on  W1E  rating  curves  (R2  increased  to  0.32)  but  removed  some  anomelously  high 
suspended  sediment  points  associated  with  moderate  discharges  on  WIN  raising  the  R2 
from  0.21  to  0.37. 
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Figure  3-19.  1999  W1  and  W1E  Sediment  Rating  Curves  for  entire 
Hydrographs  (All  Data). 
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’99  Sediment  Discharge  Rating  Curve  (w/o  1st  Flush) 


Figure  3-20.  1999  W1  and  W1E  Sediment  Rating  Curves  without  First  Flush 


Seasonal  differences  in  observed  sediment  concentrations  suggest  that  seasonal 
stratification  may  further  improve  the  rating  curves.  For  this  analysis  the  entire  1999  W1 
sediment  concentration  data  were  split  into  summer  (April  through  September)  and 
winter  (October  through  March)  samples.  The  results  of  this  stratification  improved  the 
variability  even  more  and  suggest  that  the  summer  convectional  storms  with  short 
duration  and  high  intensities  are  associated  with  high  variability  in  suspended  sediment 
concentrations  (R2=  0.38)  (Figure  3-21).  The  winter  storms  sediment  concentrations  can 
be  approximated  with  discharge  with  a  reasonable  accuracy  (R2  =  0.54).  Interpretation  of 
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these  apparent  seasonal  effects  is  complicated  by  changes  in  land-use  and  is  discussed  in 
a  later  section. 


’99  Seasonal  Sediment  Discharge  Rating  Curve 


Figure  3-21.  W1  Seasonal  Rating  Curve 
(Including  both  Rising  and  Falling  Limbs) 

Finally,  a  logical  step  was  to  see  if  removing  the  first  flush  from  the  seasonal 
rating  curves  would  yield  good  rating  curves  and  reduce  variability  to  the  seasonal 
stratification.  For  winter  storms,  however,  this  was  done  by  removing  the  data  points 
from  the  first  fifteen  minutes  following  the  precipitation  centroid.  The  results  for  summer 
storms  were  substantially  improved  from  0.38  to  0.54  explanation  of  variance  (Figure  3- 
22). 
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'99  Summer  (Convectional)  Sediment  Discharge  Rating  Curve  (w/o  1st  Flush) 
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Figure  3-22.  W1  Summer  Rating  Curve  without  First  Flush 

Not  much  improvement  for  winter  data  with  the  removal  of  the  first  flush  data  points  as 
evident  with  an  increase  of  only  6  percent  in  R2.  Winter  stratification  is  sufficient  and 
little  is  gained  by  the  first  flush  stratification.  Reduction  of  summer  variability  due  to  the 
first  flush  data  was  substantial  with  an  increase  of  R2  from  0.38  to  0.54  and  resulted  in  an 
improved  association  of  sediment  concentration  as  explained  by  discharge.  Therefore, 
seasonal  stratification  and  the  removal  of  the  first  flush  data  points  yield  an  improved  the 
summer  suspended  sediment  concentration. 

Hysteresis 

In  an  attempt  to  identify  process-related  sediment  dynamics  for  the  system,  the 
timing  of  sediment  deliveries  was  further  analyzed.  This  timing  can  be  expressed  as 
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“hysteresis.”  Hysteresis  is  the  systematic  direction  in  which  a  dependent  variable  varies 
from  a  regression  line  when  sequence  of  values  are  examined  (James,  1998).  Due  to  the 
strong  first  flush  of  sediment,  sediment  concentrations  in  WIN  typically  display  a  spiked- 
clockwise  hysteresis  (Figure  3-23).  It  is  this  spike  within  the  first  fifteen  minutes  that  is 
most  difficult  to  characterize  with  sediment  rating  curves.  Although  subsequent 
suspended  sediment  concentrations  remain  hysteric,  the  error  they  impart  is  relatively 
minor  (e.g.  points  between  26  and  31 1/s  on  Figure  3-23). 


Figure  3-23.  WIN  Typical  Clockwise  Hysteresis 


W1E  has  a  different  timing  of  suspended  sediment  concentration  than  WIN. 
W1E  does  not  have  the  quick  pulse  of  sediment  like  WIN  but  has  a  slower  response  in 
sediment  concentration  and  discharge  due  to  the  distance  of  the  channels  from  the  roads. 
Generally,  the  W  IE  sediment  concentration  peak  occurs  around  an  hour  after  the 
precipitation  centroid.  WIE’s  sediment  concentration  generally  peaks  before  or  shortly 
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after  the  discharge  peak.  This  is  illustrated  in  the  response  on  16  June  1999  where  the 
sediment  concentration  peaked  on  the  rising  limb  of  discharge  giving  a  clockwise 
hysteresis  with  the  spike  later  than  in  WIN  (Figure  3-24).  For  storms  where  W1E 
suspended  sediment  peaks  after  peak  discharge,  a  counter  clockwise  hysteresis  results 
with  the  spike  occurring  closer  to  peak  discharge.  W1E  has  varied  seasonal  response. 
With  the  quick  summer  storms,  the  sediment  concentration  response  was  about  45 


minutes  after  the  precipitation  centroid  and  the  winter  storms  had  about  a  75-minute 
delay.  This  could  be  due  to  storm  type  and  the  increased  winter  ground  litter. 


Figure  3-24.  W1E  Typical  Clock  Wise  Hysteresis 


The  response  of  W1  combines  the  varied  responses  of  W  IN  and  W  IE.  W 1  tends 
to  have  a  high  early  spike  with  the  arrival  of  the  first  flush  from  WIN  and  a  lower  peak 
with  the  later  arrival  of  sediment  from  W1E.  The  combination  of  sediment  deliveries 
from  the  two  tributaryes  generally  results  in  the  sediment  concentration  of  W1  peaking 
early  but  maintaining  a  high  concentration  and  then  having  a  smaller  second  peak  on  the 
rising  limb.  Generally,  short  duration  storms  with  the  first  flush  have  a  clockwise 
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hysteresis  and  a  very  slight  second  peak.  During  longer  events,  the  delayed  response 
from  W1E  results  in  a  smaller  but  notable  second  peak.  If  this  peak  arrives  before  peak 
discharge,  it  exaggerates  the  clockwise  hysteresis  (Figure  3-25).  If  it  arrives  after  peak 
discharge  it  causes  a  combination  of  both  a  clockwise  and  a  counter-clockwise  hysteresis 
(Figure  3-26). 


17  17.5  18  18.5  19  19.5  20 
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Figure  3-25.  W1  Clockwise  Hysteresis  with  the  Quick  WIN  First  Flush  and  the  W1E 
Sediment  Arriving  Before  Peak  Diacharge 
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W1  Single  Event  16  Jun  99 
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Figure  3-26.  W1  Complex  Hysteresis  with  Early  WIN  first  Flush  and  Late  Arrival 
of  W1E  sediment  Pulse  after  Peak  Discharge 

Unit  Flux  and  Unit  Flow  of  W1E  and  WIN 

WIN  and  W1E  rating  curves  of  unit  flux  against  unit  flow  were  created  to 
evaluate  differences  between  contributions  from  the  two  tributaries  (Figure  3-27).  By 
integrating  the  sediment  concentration  data  over  storm  durations  and  dividing  by  the 
drainage  area,  much  of  the  storm  variability  in  sediment  loads  is  removed.  For  example, 
explained  variance  is  increased  from  0.31  for  the  raw  W IE  concentration  data  to  0.90  for 
the  storm  fluxes.  A  higher  degree  of  unexplained  variability  (R  =  0.47)  in  WIN  unit 
sediment  flux  is  still  present.  ANCOVA  shows  no  significant  difference  between  W1E 
and  WIN  at  the  0.05  alpha  level  with  (F=  0.03  and  P  (>F  =  0.8574).  Therefore,  these 
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results  could  not  reject  H30  and  conclude  no  significant  difference  between  the  W1E  and 


WIN  unit  sediment  flux. 


Land-Use  Changes 

In  early  December  1999,  extensive  land-use  changes  occurred.  Active  logging 
began  with  major  portions  of  W1E  basin  being  effected.  The  channel  was  greatly  altered 
by  the  addition  of  tree  stumps  and  logs.  The  dirt  roads  were  graded  and  extensive  vehicle 
traffic  loosened  the  dirt  surface.  One  storm,  on  10  January  2000,  was  sampled  to  evaluate 
any  significant  changes  after  the  change  in  land  use.  Using  all  data  for  W IE,  reasonable 
rating  curves  were  generated  with  pre-harvest  and  post-harvest  sediment  data 
stratification  (Figure  3-28).  Unit  sediment  flux  was  different  after  harvesting  and 
ANCOVA  results  show  a  significant  difference  in  the  slope  of  the  W1E  post  harversting 
unit  flux  unit  flow  rating  curve.  WIN  had  a  uncharacteristically  high  runoff  response 
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from  this  storm.  Field  observations  indicated  that  this  atypical  response  from  WIN  was 
due  to  overbank  discharge  from  W1E  which  crossed  over  the  floodplain  about  100  m 
upstream  of  the  gage  sites  and  flowed  into  the  WIN  tributary.  This  was  the  first  time  this 
overflow  from  W1E  to  WIN  had  been  observed  and  the  high  stage  (465  mm  at  Gl)  was 
approximately  double  any  previous  observed  stage  heights  (maximum  at  215  mm  at  Gl). 
Thus  comingling  of  the  W1E  water  with  WIN  is  not  likely  to  have  effected  the  other 
samples.  This  shared  discharge  and  the  limitation  of  only  one  post-harvest  event  make  it 
difficult  to  precisely  specify  the  relative  contributions  of  runoff  for  the  two  tributaries  and 
are  viewed  as  preliminary.  The  evaluation  of  response  from  change  in  land-use  for  WIN 
was  not  done. 


Figure  3-28.  W1E  Sediment  Flux  Curves  before  and  after  Timber  Harvesting 
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SUMMARY  OF  ANALYSIS 


High  erosion  rates  and  effective  transport  of  the  sediment  on  the  dirt  roads  show 
that  the  dirt  roads  are  sources  of  accelerated  erosion  and  sediment  transport.  Much  of  the 
eroded  sediment  is  stored  at  the  base  of  the  roads  studied  but  this  may  not  always  be  the 
case.  High  variability  in  lumped  erosion  models  make  it  difficult  to  accurately  predict 
soil  loss  in  such  a  diverse  environment  especially  with  the  important  role  of  the  dirt  roads 
as  effective  channels  for  runoff  and  sediment  transport.  Alterations  to  the  soil  horizon 
through  many  years  of  road  grading  have  drastically  increased  the  difficulty  of  accurately 
parameterizing  a  model. 

Precipitation  totals  for  this  study  period  were  approximately  40  cm  below  regional 
averages  with  no  substantial  contribution  from  hurricane  events.  Stage-discharge 
relationship  at  W1E  were  established  through  field  measurements  of  ten  storm  events  and 
functionalized  with  second-order  polynomials.  At  the  W1  main  gage  the  exsiting  rating- 
curve  was  used  (Dean  et  al.  1998;  2000).  Continuous  stage  recorders  were  established  on 
both  the  east  tributary  (gage  7)  and  the  main  stream  channel  (gagel).  However, 
equipment  failures  in  May  and  December  limited  the  stage  data  recorded.  The  mV  data 
from  the  recorders  were  functionalized  to  correspond  to  field  observed  staff  gage  heights 
at  their  respective  gages  and  subsequently  converted  to  discharge  by  applying  the  stage- 
discharge  function  (Figures  2-3  and  2-4). 

Interrill  erosion  was  examined  through  2  foot  by  2  foot  wooden  small  box  plots  by 
testing  unit  sediment  mass  differences  between  dirt  road  plots  and  controls  placed  in 
forested  areas.  The  dirt  road  unit  sediment  mass  mean  was  tested  and  found  to  be 
significantly  greater  than  forested  unit  sediment  mass  at  the  0.05  alpha  level  (t-Stat  = 
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2.3152  and  P  (T<t)  =  0.0342)  (Figure  3-9).  Field  observations  indicated  that  raindrop 
splash  erosion  was  a  more  dominant  process  than  sheet  erosion  at  the  small  box-plot 
scale. 

Rill  erosion  was  measured  through  field  surveys  of  cross-sections  on  a  hillslope 
(R3,  Figure  3-7)  and  roadside  gully  (R2,  Figure  2-8).  The  R3  cross-section  five-month 
soil  loss  was  estimated  to  have  been  15.7  kg/m2  or  70  tons/ac  (Table  3-1).  The  R2  cross- 
section  five-month  soil  loss  was  estimated  to  have  been  1.42  kg/m  .  Both  erosion 
estimates  indicate  that  rill  erosion  is  dominant  in  areas  of  channelized  flow  and  resulted 
in  high  losses  of  road  sediment.  Field  observations  of  large  sediment  plumes  below  these 
sites  indicate  effective  transport  of  road  sediment  from  both  the  road  surface  and  along 
the  roadside  gullies.  Thus,  roads  produced  significantly  higher  sediment  and  Hi0  was 
rejected. 

Sediment  deposition  was  noted  at  several  locations  within  the  study  area  and 
primarily  at  locations  where  road  runoff  entered  the  forest  woodline.  Plastic  mats  were 
used  as  sediment  traps  to  measure  deposition  in  the  forest  transition  zone  or  buffer  zone 
between  sediment  sources  and  stream  channels.  The  unit  weight  (kg/m  )  of  sediment 
deposition  was  measured  across  five  transects  with  a  gradient  less  than  2  percent  (Figure 
2-9).  The  unit  weights  along  the  mat  transects  indicate  that  a  large  portion  of  the 
sediment  was  trapped  in  the  first  10  m  of  the  buffer  zone  with  a  decline  to  zero  at  around 
30  m  and  a  slight  increase  at  the  40  m  buffer  width  (Figure  3-10).  The  mean  unit  mass 
deposited  at  the  first  10  m  buffer  width  was  significantly  greater  than  the  combined  mean 
deposition  weights  at  the  other  three  sites  at  the  0.05  alpha  level  (t-Stat  =  2.27  and  P 
(T<=t)  is  0.04).  A  statistical  regression  of  the  unit  mass  against  transect  distance  explains 
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70%  of  the  variance  (Figure  3-11).  Deposition  mat  textural  analysis  reinforced  the  unit 
mass  findings  showing  that  the  greatest  percentage  of  sand  was  deposited  in  the  first  10  m 
of  the  buffer.  This  indicates  that  the  larger  sand  particles  are  settling  out  quickly  while 
the  fine  sediment  is  being  transported  further  along  the  transect.  There  is  also  an  increase 
in  organic  matter  along  the  transect  as  the  runoff  mobilizes  the  litter  layer.  Sand  textures 
measured  by  sonic-sieving  show  very  little  material  coarser  than  sand  (>  2mm)  was 
transported  and  that  the  coarse  sands  (>  0.5mm)  generally  decrease  with  an  increase  in 
buffer  width.  There  is  also  a  sharp  increase  in  the  fine  and  very  fine  sand  particles 
(<0.25mm  and  >0.63mm)  and  in  organics  (based  on  LOI)  with  the  maximum  occuring 
with  buffer  widths  of  30  m.  These  patterns  indicate  that  deposition  is  occurring  in  the 
lower  velocity  flow  environments  after  the  larger  particles  had  settled.  These  results  for 
moderate  sized  storm  event  have  significant  deposition  of  sediment  with  distance  along  a 
buffer  zone.  Thus,  H20  is  rejected. 

The  relationship  between  sediment  concentration  and  discharge  at  W1  is  complex 
due  to  the  varied  responses  from  the  two  tributaries.  The  first-flush  phenomenon 
apparent  in  WIN  is  a  result  of  road  runoff  reaching  the  stream  through  an  area  of 
inadequate  buffer  zone.  The  delayed  response  fromWIE  also  results  in  an  out-of-phase 
relationship  between  peaks  of  discharge  and  suspended  sediment  concentrations.  During 
a  high  discharge  event  on  10  Jan  2000,  the  two  tributaries  had  similar  sediment  color  and 
timing  of  peak  stages.  Field  observations  however,  found  that  this  concurrent  timing  was 
due  to  shallow  overbank  flows  passing  from  W1E  to  WIN  approximately  100  m  to  300 
m  upstream  from  the  confluence.  The  storm  sediment  flux  totals  from  five  moderate 
events  indicated  no  significant  difference  between  the  total  sediment  flux  of  W  IE  and 
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WIN  when  averaged.  Comparison  of  the  overall  unit  sediment  flux  between  the  W1E 
and  WIN  indicated  no  significant  difference  between  the  two  tributaries.  Thus,  these 
results  could  not  reject  the  third  null  hypothesis,  H30,  that  there  is  no  significant 
difference  between  unit  sediment  fluxes  at  WIN  and  W1E. 
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CHAPTER  IV 


DISCUSSION  AND  SYNTHESIS 

Results  from  this  study  show  that  soil  erosion,  sediment  transport,  and  sediment 
deposition  are  significantly  affected  by  the  presence  of  dirt  roads.  The  high  proportion  of 
dirt  road  surface  area  in  the  W 1  watershed  generated  high  rates  of  erosion,  and  the  narrow 
buffer  zone  widths  allows  these  sediments  to  be  conveyed  to  the  efficient  ditch  system. 
IMPLICATIONS  AND  FUTURE  RESEARCH 

The  presence  and  use  of  the  dirt  roads  in  the  W1  watershed,  combined  with  the 
exposed  sandy  soil  and  intense  rainfall  potential  of  the  Sandhills  environment  of  South 
Carolina,  and  inadequate  buffer  zones  combine  to  accelerate  soil  erosion  and  sediment 
transport  and  are  a  concern  for  NPS  sediment  generation.  Sediment  generation  was 
observed  along  the  road  on  the  northern  side  of  WIN  and  the  tributary  coming  from  the 
BRM-19  rifle  range  along  the  east  channel.  Soil  erosion  management  should  begin  with 
an  investigation  of  these  major  sediment  sources,  the  impact  of  dirt  roads,  and  the  ability 
of  buffer  zones  to  control  transport  of  sediment  in  this  environment. 

Erosion  and  Transport 

Erosion  from  dirt  roads  has  been  shown  to  be  larger  than  the  forested  areas;  i.e., 
H]0  is  rejected.  Although  locally  efficient  road  maintenance  and  soil  erosion  reduction 
efforts  have  been  initiated  on  Ft  Jackson,  a  relatively  small  number  of  highly  erosion- 
prone  areas  are  contributing  heavily  to  sediment  yields.  For  example,  field  observations 
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of  overland  flow  from  the  BRM-19  rifle  range  and  blocked  drainage  channels  below  it 
indicate  sustained  gully  formation  via  scouring  directly  below  this  source.  A  steel 
sediment  retention  dam  has  been  constructed  below  the  rifle  range,  but  field  observations 
and  suspended  sediment  samples  at  the  dam  outlet  indicate  substantial  sediment  loads  are 
passing  through  this  structure.  Recent  gully  activity  below  the  BRM-19  range  and 
scouring  by  overland  flow  above  the  dam  are  rapidly  producing  sediment,  which  will 
reduce  the  reservoir  capacity  and  the  useful  life  span  of  the  dam.  Reducing  overland  flow 
from  the  BRM-19  rifle  range  by  cleaning  drains  and  channels  and  stabilizing  gullies 
should  have  significant  impacts  on  sediment  production  and  sediment  yield  from  this  site. 
Field  observations  along  the  main  east  tributary  channel  revealed  that  the  majority  of 
sediment  reaching  the  channel  came  from  the  rifle  range.  Other  sources  of  sediment 
reaching  the  east  channel  were  roadside  ditches  especially  those  draining  Old  Harstville 
Guard  Road  and  Firebreak  1 1  where  it  joins  the  BRM19  rifle  range  tributary. 

Three  erosion-prone  areas  were  observed  along  the  WIN  channel.  Two  were 
areas  of  road  runoff  into  the  wetland  at  the  upstream  end  of  the  ditched  channel.  A  steep 
road  rill  that  was  shown  by  repeat  surveys  to  be  the  source  of  15.7  kg/m2  of  sediment 
during  a  six-month  period  deposited  this  sediment  at  the  base  of  the  hillslope  were  it  was 
susceptible  to  further  transport  by  subsequent  large  events.  Some  of  the  eroded  sediment 
was  carried  off-site  to  the  ephemeral  stream  network  where  it  could  be  deposited  or 
remobilized.  The  third  area  was  about  50  m  upstream  of  the  confluence  with  the  W1E 
that  drained  about  100  m  of  dirt  road  surface. 
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Deposition 


During  the  study,  flow  from  W1E  exceeded  bank-full  conditions  and  crossed  the 
floodplain  and  into  WIN,  indicating  that  in  high  runoff  events  the  wetland  loses  retention 
efficiency  and  allows  higher  sediment  yields  to  leave  the  W1E  basin  circumventing  the 
gage  site.  The  null  hypothesis  two,  that  there  would  be  no  effect  of  buffers  was  rejected. 
Thus,  buffer  zones  were  effective  in  reducing  the  amount  of  sediment  reaching  the  stream 
channel  network  during  the  moderate  magnitude  flows  monitored  by  this  study. 

Sediment  deposition  along  transects  entering  a  wetland  indicate  that  buffer  zone  widths  of 
10  m  are  effective  in  trapping  coarse  sand  and  reducing  the  transport  of  fines.  These 
results  were  only  from  the  moderate  storm  events  sampled  and  effective  widths  of  buffer 
zones  may  expand  substantially  with  larger  storms.  Areas  of  inadequate  buffer  zone 
width  below  erosion-prone  areas  may  result  in  conveyance  of  sediment  directly  to 
streams.  In  this  ditched  system,  main  channels  are  highly  efficient  at  sediment  transport 
so  these  inputs  are  quickly  delivered  to  the  system  of  ponds  downstream  of  the  W1  outlet. 
Improving  the  buffer  zone  width  by  only  a  few  meters  at  the  narrow  spots  along  the  dirt 
road  that  parallels  WIN  should  substantially  reduce  sediment  loads  in  road  runoff  that 
empties  directly  into  the  channel.  In  forest  transition  zones  below  dirt  roads  large  active 
sediment  plumes  were  observed.  These  sediment  deposits  cause  valley  aggradation  and 
may  ultimately  steepen  slopes  and  reduce  sediment  storage  capacity. 

Sediment  Yield 

High  suspended  sediment  concentrations  from  storm  events  with  high  intensities 
and  from  large  runoff  events  indicates  that  sediment  is  being  transported  from  this 
watershed.  Variability  in  sources  of  suspended  sediment  including  WIN  direct  road 
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runoff,  (the  “first-flush”  phenomenon),  and  W1E  overbank  flows  entering  WIN-make 
long-term  sediment  yield  calculations  difficult  and  somewhat  speculative  with  only  a 
small  data  set.  The  first-flush  phenomenon  was  most  apparent  in  the  north  tributary 
during  summer  (convectional  storms)  within  fifteen  minutes  of  precipitation  and  caused 
high  variability  in  sediment  rating  curves.  Removal  of  the  first  flush  decreased  the 
variability  in  sediment  discharge  rating  curves.  Seasonal  stratification  produced  rating 
curves  with  the  least  amounts  of  unexplained  variance.  Winter  storm  rating  curves  had 
nearly  the  same  variability  with  or  without  first  flush  inclusion.  Removal  of  the  first 
flush  substantially  decreased  summer  storm  sediment  discharge  rating  curve  variability. 

Sediment  loads  in  W1E  were  not  found  to  be  significantly  larger  than  sediment 
loads  in  WIN  by  analysis  of  five  storm  sediment  flux  and  runoff  discharge  values. 
Additionally,  no  significant  difference  was  found  between  unit  flux-unit  flow  rating 
curves  for  these  tributaries  through  analysis  of  covariance.  Since,  null  hypothesis  three 
could  not  be  rejected  in  both  cases  the  unit  sediment  fluxes  (scaled  per  unit  drainage  area) 
from  WIN  and  W1E  were  not  different  for  the  evens  of  this  study.  This  was  surprising 
given  the  close  proximity  of  WIN  ditch  system  and  the  greater  length  of  wetland  that 
W1E  passes  through. 

A  preliminary  geographic  information  system  (GIS)  buffer  zone  analysis  was 
conducted  to  calculate  the  dirt  road  densities  in  proximity  to  the  WIN  and  W IE  stream 
channels.  It  was  anticipated  that  the  WIN  sub-basin  would  have  a  higher  road  density 
than  the  W IE  sub-basin.  This  was  not  the  case  with  buffer  widths  less  than  fifty  meters. 
W1E  had  a  substantially  higher  density  in  the  10  m,  20  m,  and  30  m  buffer  widths  than 
WIN.  At  buffer  widths  greater  than  fifty  meters,  WIN  had  a  higher  density.  This  was 
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due  to  the  general  orientation  of  the  dirt  roads  and  firebreaks  in  relation  to  the  stream 
channels.  For  WIN,  the  firebreaks  were  mainly  perpendicular  to  the  channels  with  some 
dirt  roads  paralleling  the  channels.  W1E  had  the  opposite  orientation.  A  large 
contribution  to  W  IE  dirt  road  density  was  the  inclusion  of  streams  along  Old  Hartsville 
Guard  Road  and  the  BRM-19  rifle  range  tributary.  Additionally,  a  few  areas  of  Old 
Hartsville  Guard  Road  (approximately  100m  by  100  m)  of  timber  had  been  cleared  for  the 
timber  harvesting  machinery  and  large  trucks.  Field  observations  revealed  that  these  dirt 
roads  and  firebreaks  in  close  proximity  to  the  stream  channels  efficiently  conveyed  both 
water  and  sediment  to  the  ditched  channels. 

Future  research  is  needed  to  improve  the  suspended  sediment  data,  establish  more 
accurate  rating  curves,  and  combine  this  information  with  flow-duration  curves  to 
establish  long-term  sediment  yields  such  as  annual  loadings.  Comparisons  of  discharges 
at  Wl,  W1E,  and  WIN  indicates  that  continuity  exists  between  the  gages,  so  presumably, 
the  delivery  of  sediment  to  the  W1E  ditch  channel  by  the  BRM-19  tributary  explains  this 
observation.  Maintaining  continuous  discharge  loggers  at  the  main  channel  and  the  east 
tributary  will  allow  calculations  of  water  and  sediment  on  WIN.  These  data  will  allow 
the  development  of  future  discharge  frequency  curves,  evaluation  of  sediment  delivery 
ratios,  and  evaluation  of  the  effects  of  land-use  changes.  Improved  soil  erosion  models 
such  as  the  WEPP  model  will  enhance  the  predictions  of  erosion,  sediment,  and 
identification  of  processes  in  such  spatially  diverse  environments  (Elliot  et  al.  1994). 
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CHAPTER  V 


CONCLUSION 

The  goal  of  this  research  was  to  observe  the  physical  processes  involved  in 
hillslope  soil  erosion  and  movement  in  a  rural  Sandhills  environment.  The  bifurcation  of 
a  watershed  into  two  distinct  sub-basins  provided  the  opportunity  to  observe  their 
different  responses  in  a  paired-watershed  experiment.  Similarities  of  soils,  vegetation, 
and  climate  between  the  watersheds  allowed  study  of  the  effects  of  buffer  zones  and  dirt 
roads  on  sediment  budgets  for  the  two  watersheds.  Three  hypotheses  were  tested 
concerning  different  responses  in  erosion,  deposition,  and  transport  of  sediment.  In  short, 
the  first  two  null  hypotheses  were  rejected  but  the  third  null  hypothesis  could  not  be 
rejected. 

The  first  hypothesis,  that  entrainment  and  transport  of  sediment  along  the  dirt 
roads  is  greater  than  on  hillslopes,  was  tested  by  collecting  sediment  from  raindrop  splash 
erosion  in  small  box  plots  and  cross-sectional  measurements  along  a  dirt  road  hillslope 
profile.  Unit  sediment  masses  from  boxes  on  the  dirt  road  (49  kg/m2)  were  found  to  be 
significantly  larger  than  unit  sediment  masses  from  a  forested  area  (1  kg/m  ).  Soil  loss 
rates  were  also  determined  along  a  dirt  road  from  cross-section  repeat  surveys  and  were 
found  to  exceed  the  county  soil-loss  tolerance  factor  (T).  These  surveys  indicate  effective 
sediment  transport  along  the  roads  although  much  of  the  sediment  appears  to  have  been 
deposited  locally. 
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The  second  hypothesis,  that  a  substantial  portion  of  the  total  sediment  load  is 
stored  along  the  wetland  buffer  zone,  was  tested  by  placing  plastic  turf  mats  as  sediment 
deposition  traps  along  transects  across  buffer  zones.  Observations  of  sediment  deposition 
masses  and  particle-size  distributions  showed  that  most  deposition  occurred  within  a 
short  distance  (10  m  to  20  m)  into  the  wetland.  Sediment  deposition  along  wetland 
margins  was  significantly  greater  at  10  m  into  the  buffer  zone  than  the  combined 
sedimentation  at  the  buffer  zone  widths  from  20  m  to  40  m.  Regression  analysis  indicates 
a  sharp  reduction  occurred  in  the  transport  of  sediment  through  a  buffer  zone  that  has  a 
width  of  at  least  ten  meters  (Figure  3-1 1).  These  results  were  observed  in  moderate  size 
storm  events,  so,  sediment  transport  further  into  or  completely  across  buffer  zones  may 
occur  with  larger  storm  events.  Nevertheless,  the  efficient  detention  of  sediment  by 
buffers  during  frequent  moderate  magnitude  events  indicates  their  effectiveness  as  a  NPS 
management  tool. 

The  third  hypothesis,  that  the  north  tributary  (WIN)  will  have  a  higher  sediment 
yield  per  unit  area  than  the  east  tributary  (W1E)  was  tested  using  sediment  concentration 
data  from  water  samples  collected  during  storms.  There  was  no  significant  difference 
(alpha  =  0.05)  in  total  storm  flux  between  W1E  and  WIN  for  the  five  selected  events. 
Additionally,  ANCOVA  shows  no  significant  difference  (alpha  =  0.05)  in  unit  flux  rating 
curves  between  W1E  and  WIN.  The  lack  of  difference  in  sediment  loadings  between  the 
two  tributaries  may  be  a  reflection  of  different  responses.  During  smaller  storm  events, 
WIN  produced  higher  portions  of  unit  sediment  flux  than  did  W1E.  During  larger 
events,  W1E  produced  a  higher  unit  sediment  flux  averaged  over  the  five  storms. 
Regression  analysis  of  sediment  rating  curves  also  indicates  variable  response  in  the 
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discharge  and  suspended  sediment  concentrations  between  WIN  and  W1E.  Greater 
sediment  loads  were  expected  in  the  WIN  sub-basin  due  to  road  proximity  and  the  large 
extent  of  wetland  the  W  IE  channel  passed  through.  This  study  revealed,  however,  that 
there  is  a  very  high  road  density  in  the  lower  W1E  sub-basin  and  that  high  sediment  loads 
delivered  by  the  BRM-19  rifle  range  tributary  are  efficiently  transported  by  the  W1E 
ditched  channel. 

Recent  timber  harvesting  across  the  watershed  allowed  comparisons  of  suspended 
sediment  data  for  preharvest  and  post-harvest  periods  to  evaluate  change  in  sediment  flux 
due  to  land-use  changes.  Changes  in  land-use  were  examined  by  ANCOVA  to  evaluate 
any  response  to  changes  in  the  east  tributary.  With  sediment  flux  data  from  only  one 
post-harvest  storm  event  there  was  significant  difference  (alpha  =  0.05)  when  compared 
with  the  remaining  sediment  flux  data.  More  data  are  needed  to  determine  whether  this 
difference  is  attributed  to  seasonality,  land-use  change,  or  event  magnitude. 

In  addition  to  the  hypothesis  tests,  field  observations  and  laboratory  measurements 
revealed  several  erosion  and  sediment  processes.  For  example,  regression  analysis  of 
suspended  sediment  concentration,  discharge,  unit  flux,  and  unit  flow  determined  stage- 
discharge  functions  and  sediment  rating  curves  (Meade,  Yuzyk,  and  Day,  1990,  James 
1998).  Data  were  stratified  by  season,  type  of  precipitation  event,  rising  versus  falling 
limbs  of  hydrographs,  and  by  year  to  identify  the  dominant  processes  driving  sediment 
deliveries.  Sediment  loads  and  discharges  varied  by  tributary  and  against  previous  study 
data.  The  existence  of  a  first  flush  phenomenon,  over-bank  flooding  of  the  east  tributary 
into  the  north  tributary,  and  the  delayed  arrival  of  the  east  tributary  peak  discharge  and 
peak  sediment  loads  make  highly  variable  rating  curves.  Stratification  of  sediment 
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concentration  by  season  and  by  time  (within  the  duration  of  storms)  allowed  much 
improved  sediment  ratings  curves  to  be  achieved.  In  particular,  winter  data  alone  resulted 
in  higher  explained  variance  (R2  of  0.54).  Summer  suspended  sediment  data  (R2  =  0.38) 
had  much  unexplained  variance  and  was  further  stratified  to  isolate  the  first  flush 
sediment  arrival,  which  were  not  well  related  to  discharge.  This  increased  explained 
variance  to  (R2  =  0.54). 

This  study  indicates  that  dirt  roads  and  firebreaks  are  areas  of  accelerated  erosion 
and  act  as  part  of  an  efficient  channel  network  for  the  entrainment  and  transport  of 
sediment  particles.  The  importance  of  buffers  between  these  roads  and  stream  channels 
was  evident  in  the  significant  reduction  of  sediment  deposition  over  short  distances  across 
buffers.  Given  the  pronounced  density  of  dirt  roads  and  firebreaks  in  this  sandy 
environment,  in  areas  with  inadequate  buffer  zone  widths,  unchecked  road  runoff  acts  as 
an  important  source  of  non-point  source  pollution  to  ditches,  which  efficiently  carry  it 
down  stream. 

Future  Research 

Continued  collection  of  discharge  and  suspended  sediment  concentration  data 
from  the  east  tributary  and  the  main  channel  would  be  most  productive  in  expanding  the 
watershed  database.  Few  large  runoff  events  were  sampled  by  this  thesis  research  and 
this  would  greatly  extend  the  range  of  data  in  the  W IN  and  W IE  tributaries.  This 
expansion  of  the  data  range  would  improve  stage-discharge  functions  and  sediment-rating 
curves.  Furthermore,  the  development  of  flow-duration  curves,  from  the  continuous  stage 
data  record  would  allow  evaluation  of  the  long-term  sediment  yield.  Expansion  of 
preliminary  WEPP  model  applications  along  with  other  recent  erosional  models  would 
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allow  insight  into  the  production  and  movement  of  sediment  throughout  the  watershed. 
Historical  analysis  of  land  use  and  sedimentation  of  erosion-prone  areas  can  also  provide 
an  understanding  of  watershed  erosion  and  runoff  response  changes. 
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Appendix  A 


Unit  Sediment  Masses  from  Small  Box  Plots 


List  of  Tables 

Table  A-l  Unit  Sediment  Mass  Data 
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* 

Sample 

Mass 

Label 

(g/box) 

FI  A 

11.8 

FIB 

2.73 

F2A 

3.72 

F2D 

0.765 

F3C 

2.03 

R2B 

38.1 

R2C 

84 

R3A 

370 

R3B 

29.4 

R4A 

498 

R4B 

119 

Bold  denotes  dirt  road  sampling 
*  For  a  5  month  period  (1 6  Jun  99  to  1 9  Nov  99) 

Table  A-l  Unit  Sediment  Mass  Data 
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Road  Hillslope  and  Gully  Cross-section  Surveys 

List  of  Tables 

Tables  B-l.  R3  15m  Cross-section  Field  Data 
Tables  B-2.  R3  40m  Cross-section  Field  Data 
Tables  B-3.  R3  75m  Cross-section  Field  Data 
Tables  B-4.  R3  Length  Field  Data 
Tables  B-5.  R3  Slope  Field  Data 
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Position  length  (m) 

Height  (m) 
Date:4  Jun  99 

Height  (m) 
Date:7  Jul  99 

Height  (m) 

Date:  10  Sep  99 

Height  (m) 

Date:  19  Nov  99 

0 

0 

0 

0 

o" 

15 

0.55 

0.613 

0.753 

0.803 

35 

1.76 

1.842 

1.853 

1.9 

60 

3.61 

3.708 

3.87 

3.805 

75 

4.44 

4.508 

4.728 

4.705 

98 

5.25 

5.338 

5.658 

5.555 

90 

5 

5.188 

5.12 

Table  B-4.  R3  Length 


Position  length  (m) 

Slope 

Date:4  Jun  99 

Slope 

Date:7  Jul  99 

Slope 

Date:  10  Sep  99 

Slope 

Date:  19  Nov  99 

0 

0.0000 

0.0000 

0.0000 

0.0000 

15 

0.0367 

0.0409 

5.0263 

5.3610 

35 

0.0503 

0.0527 

5.3017 

5.4366 

60 

0.0603 

0.0619 

6.4635 

6.3545 

75 

0.0593 

0.0602 

6.3166 

6.2857 

90 

0.0556 

0.0577 

5.6981 

98 

0.0536 

0.0546 

5.7831 

5.6775 

Table  B-5.  R3  Slope 
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Appendix  C 


Depositional  Mat  Sediment  Particle  Data 


List  of  Tables 

Tables  C-l.  Depositional  Mat  Unit  Weight  Data 

Tables  C-2.  Depositional  Mat  Sediment  Particle  Analysis  Data 
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Appendix  D 


Suspended  Sediment  Concentration  and  Discharge  Data 
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Table  D-1.  12  Feb  1999  Single  Storm  Data 
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Table  D-2.  16  Jun  1999  Continued  Single  Storm  Data 
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Table  D-4.  19  Jul  1999  Single  Storm  Data 
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Table  D-8.  4  Oct  1999  Single  Storm  Data 
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Date  Time  Elapsed  Sample  Suspen  Sed  Stage  (6)  Flow  Suspended  Elapsed  Unit  Unit  Suspended 

Time  Concentration  Sediment  Flux  Time  Flow  Sediment  Flux 

(dec  hrs)  (mg/I)  (mm)  (L/s)  (kg/d)  (min)  (L/s/km2)  (kg/d/km2) 
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25-Jun-99  915  0.60  3  67.67  115  6.65  38.89  36  6.10  35.68 

25-Jun-99  925  0.77  4 _ 20.83 _ 1U _ 6J32 _ 11-92 _ 46  6.08 _ 10-94 

Table  10.  W1  N  Unit  Flux,  Unit  Flow,  and  Hydrograph  Data 
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Date  Time  Elapsed  Sample  Flow  Suspen  Sed  Stage  (7)  Elapsed  Suspended  Unit  Unit  Suspended 
Time  Concentration  Time  Sediment  Flux  Flow  Sediment  Flux 

_ (dec  hrs) _ (L/s) _ (mg/I) _ (mm) _ (min) _ (kg/d)  (L/s/km2)  (kg/d/km2) 

W1E  12  Feb  99 
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25-Jun-99  904  0.42  2  8.87  41.24  297  25  31.59  1.61  5.73 

25-Jun-99  915  0.60  3  9.03  40.00  298  36  31.19  1.64  5.66 

25-Jun-99  925  0.77  4  9^5 _ 40.00 _ 300 _ 46 _ 32.31  1.70 _ 5J36 
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16-Jun-99  1820  2.10  17  126  40.91  94.74  196  334.84  6.20  50.73 

16-Jun-99  1827  2.22  18  133  39.40  74.47  193  253.49  5.97  38.41 

16-Jun-99  1834  2.33  19  140 _ 38.91 _ 72.29 _ 192 _ 243.00  5.89 _ 36.82 
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Table  12  Continued.  W1  Unit  Flux,  Unit  Flow,  and  Hydrograph  Data 
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SS  Flux  kg/d 
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19.50 


Figure  D-1 .  Single  Storm  Event 
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Figure  D-3.  Single  Storm  Event 


19  Jul99 


Figure  D-4.  Single  Storm  Event 
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Figure  D-5.  Single  Storm  Event 


Figure  D-6.  Single  Storm  Event 


Figure  D-7.  Single  Storm  Event 
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SS  Flux  (kg/d) 
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Figure  D-8.  Single  Storm  Event 


10  Jan  2000 


Elapsed  Time  (min) 


Figure  D-9.  Single  Storm  Event 
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Continuous  Data  Logger  Discharge  Data 
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E  Aug  to  Dec  1 999  Hydrograph 
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